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ABSTRACT 

Detection of gravitational waves from accreting neutron stars (NSs) in our galaxy, due 
to ellipticity or internal oscillation, would be a breakthrough in our understanding 
of compact objects and explain the absence of NSs rotating near the break-up limit. 
Direct detection, however, poses a formidable challenge. Using the current data avail- 
able on the properties of the accreting NSs in Low Mass X-Ray Binaries (LMXBs), 
we quantify the detectability for the known accreting NSs, considering various emis- 
sion scenarios and taking into account the negative impact of parameter uncertainty 
on the data analysis process. Only a few of the persistently bright NSs accreting at 
rates near the Eddington limit are detectable by Advanced LIGO if they are emitting 
gravitational waves at a rate matching the torque from accretion. A larger fraction 
of the known population is detectable if the spin and orbital parameters are known 
in advance, especially with the narrow-band Advanced LIGO. We identify the most 
promising targets, and list specific actions that would lead to significant improvements 
in detection probability. These include astronomical observations (especially for un- 
known orbital periods), improvements in data analysis algorithms and capabilities, 
and further detector development. 

Key words: accretion, accretion disks — gravitational waves — stars: neutron — stars: 
rotation — X-rays: binaries — X-rays: bursts 



1 INTRODUCTION 

A number of interferometric gravitational wave observato- 
ries have been built with the intention of opening a new ob- 
servational window for studying astrophysical objects. These 
include the LIGC0, GEC0, VIRGC0 and TAMA0 detec- 
tors. The TAMA and LIGO Scientific Collaborations have 
demonstrated their ability to reach sensitivity goals, take 
year-long stretches of data with good duty cycles, and anal- 
yse the data to set, with confidence, upper limits on the 
emission from a number of possible sources. At present the 
large LIGO and VIRGO detectors are performing a signif- 
icant sensitivity upgrade. This is therefore a good time to 
carry out a realistic study of the challenges of searching for 
gravitational waves from one class of sources that future up- 
grades may render detectable: spinning neutron stars (NSs) 
in Low Mass X-ray Binaries (LMXBs). 



1 http://www.ligo.caltech.edu 

2 http://gco600.aei.mpg.de/ 

3 http:/ /wwweascina. virgo.infn.it 

4 http://tamago.mtk.nao.ac.jp 



For the la st ten years, ever since a paper by one of us 
(|Bildstenlli998h suggested that LMXBs could be steady bea- 
cons of gravitational waves, the study of LMXBs (and espe- 
cially of the source Scorpius X-l) has been one of the scien- 
tific goals of the development of detectors with greatly im- 
proved sensitivity. Since then, X-ray astronomers have gath- 
ered a wealth of new data on these sources, astrophysicists 
have built better models, and gravitational wave scientists 
have gained considerable experience of their ability to ex- 
tract weak signals from data streams. It is timely, therefore, 
to revisit these estimates. 

The expectation that LMXBs could be strong steady 
sources of gravitational waves originates in one of the 
most important outstanding questions about NSs: why 
they all spin at frequencies much less than the break-up 
limit of 1 kHz |Lattimer "fc Prakash 200^). Simple estimates 
()Cook. Shapiro fc Teukol skv 1994) of the spin-up timescales 
for accreting NSs in LMXBs suggest that there should be no 
difficulty in reaching at leas t 1 kHz. The fastest r otating ac- 
creting pulsar is at 599 Hz (iGallowav et al.ll2005 ) and many 
NS rotate much more slowly (jChakrabartv et alj|2003 ). The 
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millisecond radio pulsars, the likely offspring of LMXBs, 
also spin at slowe r rates than expect ed. The current record 
stands at 716 Hz ()Hessels et al.ll200ri ) and surveys in nearby 
glob ular clusters continue to reveal a paucity of rapid rota- 
tors jRansom et al.j|2005l ; iHessels et al] 120071 ; iFreire et al~l 
120081 ). Whether the apparent spin limit is genuine remains 
to be resolved. Nevertheless, given the current sample, it 
certainly appears that there is some brake that prevents ac- 
creting NSs from reaching the break-up limit. 

The candidate mechanisms fall into two main camps. In 
the first, accretion torques are reduced and eventually bal- 
anced by the interact ion between the accre t ion disk and the 
NS's magnetic field dchosh fc Lamblll978l : Iwliite fc Zhang] 
1 19971 : lAndersson et all 20051 ). The second possibility, and 
our focus here, is the loss of ang ular momentum via the emis- 
sion of gravitational radiation (jPapaloizou fc Pringldfl97l 
IWagoner|[l98i lBildsterJ[l99sh . There are many conceivable 
ways for an accreting NS to develop a quadrupolar asymme- 
try that lea ds to gravitatio n al wave (GW) emission: crustal 
moun t ains dBildstenl 19981 : Ushomirskv. Cutler fc Bildstenl 
200d: iMelatos fc Pavnd 120051: iPavne fc Melatosl 1200a 
Haskell. Jones fc Anderssonl l2006i ) , magnetic deformation 
ilCutlerj|2002i; lHaskell et all 120081) or internal r-mode oscil- 
lations dBildstenl 1 19981: lAndersson. Kok kotas & Stergioulasl 
iLevinl 1 19991 : lAndersson et al.1 l200d: 



2002 



lAndersson. Jones fc Kokkotasl 
20021: iNavvar fc Owenl 120061 ; 
Bondarescu. Teukolskv fc Wasserman 20071 ). The exciting 
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prediction from all GW emission scenarios is the possible 
direct detection of an accreting NS by a ground-based 
interferometric gravitational wave detector. 

Searches for periodic GWs from NSs have already been 
performed with the LIGO and G EO detectors. These include 
search es of known radio pulsars (jAbbott et al.ll2004l . l2005bl . 
l2007bl ). assuming the phase of the GW signal to be locked to 
the known NS rotat ion, and wide parameter surveys f or hith- 
erto unknown NSs (jAbbott et al.ll2005al . l2"007al . l2008l ). There 
have also been search es for GWs from the a ccreting NS Scor- 
pius X-l. The first (jAbbott et al.1 l2007al ) used a coherent 
statistic on a template grid utilising 6 hours of data (limited 
by computational requirements) from the second science run 
of the two LIGO 4 km interferometers, followed by a coinci- 
dence analysis bet ween candidates from the two detectors. 
The second search (jAbbott et al.ll2007d ) used a method of 
cross-correlating the outputs of the two LIGO 4 km detec- 
tors using about 20 days of coincident data taken during the 
fourth science run. Data from resonan t bar detector s have 
also b een used in these searches; see e.g. lAstone et all (|2002l . 
120051 ) for a blind all-sky search in a narrow frequency band 
using data from the EXPLORER detector. 

These early searches, while so far yielding only upper 
limits, have served to develop and prove data analysis meth- 
ods that will be used on data from future, more sensitive 
searches. The VIRGO detector is approaching its design sen- 
sitivity and the LIGO detectors have just completed a full 
run at their first-stage design sensitivity. An initial upgrade 
of LIGO is about to commence, followed by a full upgrade 
to Advanced LIGC0 in the next decade. VIRGO expects to 
parallel these developments. These impending improvements 



http: / /www.ligo.caltech.cdu/advLIGO 



in GW sensitivity, coupled with the experience of perform- 
ing realistic data analysis, make the time right to assess 
what the relevant accreting NS properties tell us about the 
ultimate detectability of NS sources. 

This paper explains the prospects for detecting GWs 
from accreting NSs, and identifies the most promising tar- 
gets. We start in Section [2] with the best case scenario, 
highlighting the data analysis challenges and the need for 
detailed information on the accreting NS properties, which 
are presented in Section [3] Section [4] explains fully the GW 
search data analysis challenges and the direct consequences 
of parameter uncertainty on integration times and detection 
statistics, allowing us to assess future detectability in Sec- 
tion [5] We close in Section [6] by highlighting where progress 
can be made in the short-term on NS source properties, and 
the implications for current and future gravitational wave 
searches. 



2 OVERVIEW AND BEST CASE ESTIMATES 

The potential strength of a periodic GW signal at frequency 
v from s pecific LMXBs with accretion rates M was esti- 
mated by iBildstenl (|l998l ) (hereafter B98) in the mountain 
scenario. In this model a quadrupole moment Q, that is sta- 
tionary in the rotating frame of the star (a 'mountain'), leads 
to the emission of gravitational waves with a predominant 
frequency v = 2ia,, v s being the spin frequency of the NS. 
The assumption is that the accretion torque on a NS of mass 
M and radius R 



N a = M (GMR) 



1/2 



is balanced by the GW torque 

32GQ 2 ^ 



N gw = 



5c 5 



(1) 



(2) 



where fl s — 2nv s is the angular frequency of the star. The 
accretion rate is estimated from the bolometric X-ray flux F 
by assuming that the luminosity is L m GMM / R, yielding 



M = 



AnRd 2 F 
GM ' 



(3) 



where d is the distance to the source. The pre- 
dicted gravitational wave amplitu de ho (as defined in 
Ijaranowski. Krolak fc Schutzl (|l998l )1^l can be written in 
terms of the GW luminosity, E gv , 

5G 



, 2 
h 



Under the condition of torque balance, 



(4) 



h = 3x10 F_ 



27 rpl/2 



( R 



VlOkm 



3/4 /l.4A/ e V /4 /lkHz^ 1 ' 
M 



(5) 

where F-g — F/10 8 erg cm 2 s 1 . Equation ([5]) makes it 
clear that the GW signal strength depends on two observ- 
ables , the flux on the sky from the LMXB and the NS spin 
rate (jWagonerll 19841 : iBildstenll 19981 ) . This amplitude is then 



6 Note that flux formulae are often given in terms of an angu- 
lar and time-av e raged ampli tude h rather than in terms of hp. 
In lOwen et al.l dl998l) and lAndersson. Kokkotas h Stergioulasl 
for example, the quoted amplitude h? = 2/iq/5. 
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compared directly to the best case detectable amplitude in 
a long search, h^ ns , which we will derive shortly. 

These amplitudes are sufficiently weak that long 
stretches of data must be folded, using predicted signal tem- 
plates. This is called matched filtering, and its sensitivity 
improves with the square-root of the observation time, pro- 
vided that the template manages to keep phase with the real 
signal to within about one radian over the entire duration. 
Where spin and orbital parameters are poorly constrained, 
many templates must therefore be searched. Analysis can 
become statistically and computationally untenable, and 
computational loads may limit integration times and thereby 
detectability. Contrast for example the most recent searches 
for ra dio pulsars, where a p recise timing ephemeris is avail- 
able (|Abbott et al.l l2007bT ). with se arches for the closes t 
accreting neutron star Scorpius X-l (| Abbott et al.l [2007c! ) . 
Poor constraints on the spin and orbital parameters for Sco 
X-l necessitate multiple templates. This reduces the feasible 
integration time substantially, resulting in upper limits that 
are approximately two orders of magnitude larger than those 
obtained for the best radio pulsar. We begin our discussion 
there. 

2.1 Spin Frequencies of Accreting NSs 

The discoveries and studies with the Rossi X-ray Timing 
Explorer (RXTE) over the last ten years have dramati- 
cally improved our understanding of stellar spin rate, v s 
jStrohmaver fc Bildstenl 120061 : Ivan der Klisl l2006f ). Having 
some constraint on the spin is, as will become clear in later 
Sections, the most important factor determining the feasibil- 
ity of gravitational wave searches. For this reason we concen- 
trate only on sources for which there is some measurement 
or estimate of v s . We now summarize the three relevant cat- 
egories of accreting NSs: accreting millisecond pulsars, burst 
oscillation sources, and kiloHertz Quasi-Periodic Oscillation 
(kHz QPO) sources. 

The spin is measured directly in the accreting millisec- 
ond pulsars, where fixed hotspots, presumably at the mag- 
netic footpoints, are a permanent asymmetry. In 1998 only 
one such object was known; there are now ten members 
of this class, three of which show only intermittent pulsa- 
tions. The other measures of spin are indirect. Probably the 
most reliable are burst oscillations, seen during Type I X-ray 
bursts (when accreted material burns in an unstable ther- 
monuclear flash). In 1998 there were six burst oscillation 
sources: there are now twelve stars with burst oscillations 
seen in multiple bursts (including three pulsars), and seven 
stars with tentative detections. 

For the three accreting millisecond pulsars that 
also show burst oscillations, the burst oscillation fre- 
quency is at or very close t o the known spin fre- 
quency dChakrabartv et al.1 120031 : IStrohmaver et all 120031 : 
lAltamirano et alj I2008T ). It would therefore seem reason- 
able to equate the burst oscillation frequency with the 
spin frequency for the non-pulsing sources (as was done 
in B98). The frequency for a given source is highly con- 
sistent from burst to burst, implyi ng that there is at least a 
strong dependence on stellar spin (|Strohmaver et al"1ll998al : 
iMuno. Ozel fc Chakrabartvll2002l ). The detection of highly 
coherent oscillations lasting several hundred seconds dur- 
ing a superburst adds further support to this hypothesis 



jStrohmaver fc Markwardtll2002i ) . There are, however, some 
complicating factors. Firstly, burst oscillations can exhibit 
frequency drifts of up to a few Hz. Secondly, there are some 
differences in the properties of the burst oscillations of the 
non-puls ing LMXBs as comp a red t o the pulsars (see for 
example IWatts fc Strohmaverl (|2006l )). suggesting that the 
mechanism may differ. In models that involve global modes 
of the surface layers, for example, the observed freq uency 
would be offset from spin frequency by several Hz (|Hevll 
12004 IPiro fc Bildstenll2005l ) . 

The third class of sources to be considered are those 
which exhibit twin kHz QPOs. Early observations suggested 
that although the frequencies could shift, their separation 
remained relatively constant, implying a link to stellar spin. 
We now know that separation varies (often quite substan- 
tially) as accretion rate changes, and the cause of the kHz 
QPOs is still not understood. However, most models still 
depend in some way on the stellar spin, either directly or 
via the influence on t he space-time in the inner regions of 
the accretion disk (see Ivan der Klisl (|2006l ) for a recent and 
comprehensive review). Observations of twin kHz QPOs in 
two of the accreting millisecond pulsars have done little to 
resolve the situation: in one case separation is rather close to 
the s pin frequency; in the other it is slightly less th an half the 
spin (|Wiinands et alj|2003al : iLinares et al.ll2005l) . Whether 
there is any link between kHz QPO separation and stellar 
spin will doubtless emerge in due course. However, given 
that most models still predict some relation, we follow B98 
and include these stars in our analysis. 



2.2 Optimal Gravitational Wave Detection 

Before being scaled down for their current sensitivity up- 
grades, the first generation LIGO interferometers operated 
for nearly two years at or better than their design sensitivity, 
by the end of which the first generation Virgo detector was 
not far behind. The "Enhanced" LIGO and VIRGO interfer- 
ometers should begin operation in 2009, and the Advanced 
LIGO interferometers are now funded and should be opera- 
tional by 2014. Advanced VIRGO is expected on the same 
schedule. Beyond this, there are ambitious plans for a third 
generation Einstein Telescope (ET) in Europe. 

The sensitivity of a gravitational wave detector is de- 
termined by the power spectral density of its instrumental 
strain noise, normalized to an equivalent gravitational-wave 
amplitude h(t). This is just the Fourier power spectrum of 
h(t) and is called Sh(y)- the noise power per unit frequency. 
It is conventional to plot h{v) = [Sh^v)] 1 ^ 2 , which allows 
comparison of the noise to the signal's amplitude if one 
knows the signal's bandwidth. The sensitivity curves of cur- 
rent and future detectors are shown in Figure [1] The design 
noise curves for Initial LIGO, Virgo and ET are taken from 
analytic models available, for example, in the LIGO software 
repositorjQ. The most recent science runs of the LIGO de- 
tectors have in fact reached their design goals over a broad 
frequency range of interest, above ~ 40 Hz. The Enhanced 
LIGO noise curve in Figure [T] is a realistic estimate of what 
can be achieved; thus at low frequencies (below ~ 40 Hz), 



http : //www. lsc-group .phys .uwm.edu/daswg/projects/lal . html 
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the Enhanced LIGO noise curve lies above the Initial LIGO 
design curve. 

The Advanced LIGO detector configurations have not 
yet been finalized, and we therefore need to consider differ- 
ent possibilities. Advanced LIGO can potentially be oper- 
ated in a narrow band mode where sensitivity is gained in 
a relatively narrow freque ncy range at t he expense of broad 
band sensitivity (see e.g. iMeersI (|l98cf ); iBuonanno fc Chenl 
This could be particularly relevant for periodic sig- 
nals where the frequency is well known. The first panel of 
Figure [1] shows an example of such a narrow band noise 
curvfjf] The noise curve can be tuned to target a broad range 
of frequencies by changing a number of interferometer pa- 
rameters, and for our purposes, the value of the noise at the 
minimum is especially important. Thus, the lower panel of 
Fig. \T\ shows the lower envelope of the narrow band noise 
curves above 100 Hz. Between 100 Hz and about 400 Hz, 
the narrow band curves are limited by thermal noise, and by 
quantum noise at higher frequencies. See Section fo.-ll for ad- 
ditional discussion. It is important to keep in mind that the 
narrow-band "envelope" does not represent any particular 
interferometer configuration but is rather a superposition of 
many configurations. It is only useful for targeting narrow- 
band signals with frequency uncertainties of, say, O(10) Hz 
in which case we can choose the appropriate element from 
the set of configurations used to produce the envelope. Fi- 
nally, note that possible designs for the third generation de- 
tectors are still being explored, and thus the noise curve for 
the Einstein Telescope is much more preliminary. 

To start addressing the question of detecting GWs from 
LMXBs, we start by asking how strong the signal would 
need to be for detection if we knew all source parameters 
to sufficient accuracy that only one template was needed 
(a coherent fold). Computational cost is no issue, and we 
can easily integrate for long periods T b s . Such a fully co- 
herent search with D detectors of comparable sensitivity 
is a best case and defines a signal-to-noise ratio (SNR) as 
SNR 2 = hlT ohs D / S n ; thus, the SNR-squared builds up lin- 
early with the observation time (here S„(v) is the power- 
spectral density of the detector noise). Conversely, for a 
given choice of SNR threshold for detectability, this leads 
to a minim um detectable signal amplitude hp (f ollowing the 
notation of Ijaranowski. Krolak fc Schutj (|l998l )): 



ho 



11.4 



S n. 



(6) 



DT hs 

which is a useful indicator of the sensitivity of such a search 



8 The various Advanced LIGO noise curves shown here have 
been obtained using version 6.2 of the Matlab script "Bench". 
It is worth mentioning that a newer version of Bench, v7.0, is 
currently under development which includes an improvement in 
coating thermal noise; it is probably worth re-calculating the Ad- 
vanced LIGO noise curves once this and later versions become 
available. 

9 To generate the lower envelope of the narrow band noise curves, 
we have used Bench to calculate the narrow band noise curves for 
a range of choices for a few interferometer parameters. These pa- 
rameters are the phase and transmittance of the signal recycling 
cavity, and the transmittance the input test mirror and the power 
recycling mirror. The envelope is then obtained by calculating the 
convex hull of the minima of the various noise curves. 
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Figure 1. The noise curves for Initial LIGO, Virgo, Advanced 
LIGO, and the third generation Einstein Telescope (ET) interfer- 
ometers. The Initial, Enhanced and Advanced LIGO curves are 
labelled "I-LIGO", "E-LIGO" and "A- LIGO" respectively; "ET" 
is Einstein Telescope. There arc three curves for Advanced LIGO 
in each panel. The top panel shows a nominal broad band config- 
uration (the so called "zero de-tuned" configuration) , an example 
of a narrow band curve ("A- LIGO NB EX"), and the total ther- 
mal noise ("A-LIGO TH"), i.e. the sum of the suspension and 
mirror thermal noise curves. The lower panel shows the lower en- 
velope of the narrow band curves ( "A-LIGO NB ENV" ) instead 
of the narrow band example. The thermal noise is shown because 
it is sometimes taken as a theoretical lower bound on the narrow 
banding for frequencies above, say, 100 Hz; as seen from the lower 
panel, this is not a good approximation at higher frequencies. 



The factor of 11.4 corresponds to an SNR threshold which 
would lead to a single tria l false alarm rate o f 1% and a 
false dismissal rate of 10% (jAbbott et al.ll2007al ) and a uni- 
form aver aging over all possible source orient ations and sky 
positons (|jaranowski. Krolak fc Schutzlll998h . 

Figure [2] compares the predicted and detectable am- 
plitudes for this best-case scenario using the long-term flux 
average derived in Section [3] and summarised in Table [TJ and 
assuming that each NS is in perfect spin balance (so that we 
can neglect spin derivatives) , with gravitational wave torque 
balancing that of accretion. We assume that the spin fre- 
quency for the kHz QPO sources lies in the middle of the 
known range of separations, and take T ^ s = 2 years as a 
reference value. 

Figure [2] naively implies that the kHz sources are the 
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Figure 2. Best case detectability for the mountain scenario for T b s = 2 years, balancing the long-term average flux, and assuming 
that all parameters are known (single template search). The bursters are divided into two groups: those for which the frequency is 
confirmed (filled) and those for which the frequency requires confirmation (open), see Section [3] The frequency at which the kHz QPO 
symbols appear is derived from the centre of the measured range of separations: the predicted amplitude would be higher if the frequency 
were lower, and vice versa. We show detectability thresholds for Initial LIGO (I-LIGO), Enhanced LIGO (E-LIGO), Advanced LIGO 
(A-LIGO), and the Einstein Telescope (ET). We also show two detectability curves for Advanced LIGO Narrow Band: the expected 
envelope for the narrow band detector that includes all sources of noise (A-LIGO NB), and a curve showing only the thermal noise floor 
(A-LIGO NB TH). 



most easily detected sources. However, they offer a specific 
challenge, as we do not know many of their spins or orbital 
parameters. Detecting the GW signal requires a knowledge 
of the GW phase evolution which depends crucially on the 
orbital parameters. Ignorance or inaccurate knowledge of the 
orbit could then require a search over a significant number 
of parameters and as we shall see, this can have a dramatic 
impact on the sensitivity of the GW search. Hence assessing 
the detectability of GWs requires a more careful descrip- 
tion of the mechanics of the gravitational wave data analysis 
process. Most of this paper is therefore focused on clearly 
assessing the data analysis challenges for present and future 
generations of GW detectors that takes into account the 
limitations imposed by incomplete astrophysical knowledge 
and finite computational resources. 



3 ACCRETING NEUTRON STAR 
PROPERTIES 

In this Section we summarize our knowledge of the rel- 
evant search parameters for all accreting NSs in LMXBs 
where there is some estimate of the spin: the accret- 



ing millisecond pulsars, the burst oscillation sources, and 
the kHz QPO sources. This sample contains both atoll 
and Z sources, a classifica tion determined by the spec- 
tral and timing properties (lHasinger fc van der Klisl Il989l; 



Muno Remillard fc Chakrabartvl |2002| ; iGierlinski fc Pond 
2002 : Ivan der Klisl l200q ) . As explained in Section [2] we re- 
quire flux histories, spin frequencies and orbital parameters, 
some of which are measured, whilst others can only be esti- 
mated. We also need to gauge the uncertainty on each quan- 
tity. This sets the required search parameter space, and in 
Section[5]we use this information to compute the number of 
signal templates required for each source. 

This Section is rather lengthy, as we give full details of 
the provenance of all of the values used in our study. There 
are two main reasons for this: firstly, to make clear the link 
between the astronomical observations and the consequences 
for gravitational wave searches. Secondly, many of the val- 
ues that we derive involve assumptions, or draw on old or 
uncertain measurements: where this is the case we wanted 
to make it explicit, in order to drive future astrophysical 
modelling and observations. The information in this Section 
should also be a useful resource for anyone intending to carry 
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out a search for gravitational waves from these objects, or for 
future detectability studies. Readers who are not concerned 
with the details of the source properties and uncertainties 
can however skip this Section, and will find the key results 
summarized in Tables [T] - [4] 



3.1 Constraining Fluxes and Accretion Rates 

The observed flux, F, sets the expected GW signal strength 
(Equation[5}. For persistent sources we record only the long- 
term average flux, F av , whereas for transient sources we also 
record the outburst flux, F b (averaged over outbursts and 
quiescent periods). Bolometric outburst fluences have been 
computed for the accreting millisecond pulsars, but for the 
majority of sources this is not the case. For the remaining 
sources we determine F using data from the RXTE All Sky 
Monitor (ASM), which provides a near-continuous history of 
source activity from 1996 to the present in the 2-10 keV band 
jLevine et al.ll 19961 ). Converting the 2-10 keV flux to a bolo- 
metric flux requires detailed spectral mod elling. For most 
burst oscillation sources we use the results of lGallowav et all 
|2008l ): these authors carry out spectral modelling using the 
pointed RXTE PCA data (2.5-25 keV) to estimate bolo- 
metric flux. We compare the calculated fluxes to the ASM 
countrate at the time of the observation to establish this 
relationship. For the transient sources we only include in 
our integrated ASM histories the times when the source is 
in outburst and detectable above a 3<r level. This avoids 
contamination from other sources in the field which would 
overestimate the long-term flux. This method of estimating 
flux history introduces some errors, since spectral shapes 
(and the correction from 2-10 keV flux to bolometric flux) 
will change: but it gives a reasonable estimate. This uncer- 
tainty should however be borne in mind in Section [6] for 
those sources that are on the margins of detectability. 

We also record position and distance. Source position 
must be known to a certain degree of precision for the long 
folds that this type of analysis necessitates ( Section I4.ip . We 
have listed the most accurate and up to date position known 
for the X-ray source or its optical, infra-red or radio coun- 
terpart. Source distance, which is relevant for the emission 
modelling in Section [5] can be estimated in several different 
ways. Only for the closest source in our sample, Sco X-l, 
can the distance be measured via parallax. For the other 
sources different methods are used: location of the source 
in a globular cluster; the presence of radius expansion X-ray 
bursts (where luminosity reaches the Eddington limilF°h: in- 
ferences about mass transfer from the long-term X-ray flux, 
assuming that the binary orbit evolves due to gravitational 
wave emission; absorption and spectral modelling. 

We give details of pulsar frequency, burst oscillation 
frequency, and twin kHz QPO separation, as measured with 
RXTE's Proportional Counter Array (PCA). For sources 
with only a burst oscillation frequency, we assume that v a 
lies within ±5 Hz of the burst oscillation frequency Vb- For 
the kHz QPO sources we assume that the spin lies within 



10 The Eddington limit depends on the composition. Sometimes 
composition can be inferred from burst properties, but this is 
not always the case, leading to additional uncertainty. See also 
iGallowav. Ozel fc Psaltis] 1120081 1. 



the reported range of kHz QPO separations. We review both 
of these assumptions in more detail in later Sections. 

3.2 Orbital Uncertainties 

We must know (or presume) the orbital parameters in or- 
der to fold long stretches of gravitational wave data, most 
importantly, the orbital period P rb and eccentricity, e. We 
also need a reference time within the orbit. Depending on 
the original reference we quote either Tgo, the time of infe- 
rior conjunction of the companion star, or T asc , the time of 
ascending node, when the Doppler shifted frequency of the 
neutron star is at its lowest. Note that Tgo = T asc + P or b/4. 
The third parameter is the projected semi-major axis of the 
neutron star, a x sini, which we denote a p . Depending on 
the measurement, we may instead quote the amplitude of 
the projected orbital velocity of the neutron star, « x sini 
(referred to as Kl in the optical literature). The two quanti- 
ties are related by v x sini = 27ra x sin i/P or b • The orbital pa- 
rameters are measured directly for the accreting millisecond 
pulsars via X-ray timing. The situation is more challenging 
for the non-pulsing and intermittent sources. 

Some high inclination systems show eclipses in the X- 
ray lightcurve, providing both P or b and T asc . There are also 
systems that show dips rather than full eclipses: the dips 
occur when the NS is obscured by the bulge where the 
mass stream from the donor star joins the accretion disk 
(IWhite fc Swan k1 ll982T ). Dips certainly tell us P or b, and re- 
strict the inclination to lie in the range 60° — 75° . What they 
do not necessarily yield is T asc : in the two systems that show 
both dips and eclipses, dips occur at various offsets from the 
eclipse times dCominskv fc Wood! 19841 : iParmar et al.|[l986l : 
iMotch et all Il987l ; ISmale et al.l I1992T ). Detailed modelling 
is therefore required to determine the relationship between 
T asc and the dip time Tdi p . Throughout this Section we list 
the most recent orbital ephemeris for each system. In sev- 
eral cases the ephemerides are sufficiently out of date that we 
should consider T asc to be totally unconstrained. However, 
we presume that it would be straightforward to obtain a new 
ephemeris, with an error no worse that that of the existing 
measurement. We therefore use all existing measurements in 
our initial assessment of detectability. 

The orbital parameters can also be measured in wave- 
bands other than the X-ray, particularly the optical. A num- 
ber of systems show photometric variability at the orbital 
period. Maximum optical light occurs when the NS is at 
inferior conjunction and we observe reprocessed emission 
from the heated face of the donor star, yielding both P or b 
and T asc , but not the projected semi-major axis. An alter- 
native method that can provide all three pieces of infor- 
mation is phase-resolved optical spectroscopy. LMXBs ex- 
hibit many emission lines, some fro m heated face of the 
donor star jSteeghs fc Casaresl [20021 ) . others from the ac- 
cretion disk very close to the compact object. By mea- 
suring the orbital Doppler shifts of these lines, and using 
techniques such as Doppler tomography to check the emis- 
sion location, major progress has been made in comput- 
ing orbital parameters. Additional constraints on the sys- 
tem are possible (assuming Roche lobe overflow and tidal 
locking) if one can detect rot ational broadening of absorp- 
tion lines from the donor star (IHorne. Wade fc Szkodvlll986l : 



ICasares. Charles fc Kuulkersll 19981) . 
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Unfortunately, for many of the systems of interest, one 
or all of the orbital parameters are unknown. We can how- 
ever still place various constraints on the systems to reduce 
the number of templates required. For those systems with 
measured orbital periods, for example, we can assume that 
the donor star fills its Roche lobe. This fixes pa, the mean 
density of the donor star: 



(7) 



The range of donor types with this p limits the donor mass 
Ma = m^MQ. LMXB donors include main sequence stars, 
evolved hydrogen-burni ng stars, helium-burning stars, and 
brow n or white dwarfs (|Podsiadlowski. Rappaport fc Pfahll 
120021 ). For Porb > 10 hours, pa < 1 g cm 3 , the companion 
must be an evolved hydrogen-burning star such as a sub- 
giant. For ultracompact systems with P orD < 80 minutes, 
the companion must be a white dwarf, a helium star, or 
a highly evolved helium-rich secondary dRappaport fc Jossl 



1984 iDelove fc Bildsten|[200l ; iNelemans. Jonker fc Steeghsl 



20061 ). For intermediate orbital periods more options are pos- 



sible, maximum donor mass being by the main sequence star 
- but evolved, les s massive companions or even brown dwarfs 
may be possible dTout et al.ll 19961; IChabrier fc BaraffelpOOOl ; 



IChabrier et~aTl |2000| ; iBildsten fc ChakrabartvibOOlH 11 !. Ad- 
ditional constraints on donor properties may come from 
X-ray burst properties or spectral type: more evolved, 
lower mass, donors will have later spec tral types than 
the main sequence star with the s ame pa (jBaraffe fc Kolbl 
l200d ; iKolb. King fc Baraffd 1200 ll ). Having finally estab- 
lished the range of likely donor masses, and knowing that 
the NS mass M x = m x M B is in the range 1.2 — 2.4M© 
(ILattimer fc Prakashl [20071 ). we then estimate a x sin i using 



a x sini = 1.174 



Porb 

Th7 



2/3 



(m d + ra x ) 2 / 3 



sin i It — s . (8) 



For systems that do not show dips or eclipses (i < 60°), 
and where there is no other limit on inclination, Equation 
([8| gives only an upper limit. For dipping systems we as- 
sume 60° < i < 75°, giving both upper and lower limits. 
For eclipsing systems an additional c onstraint comes from 
the h alf-angle of the X-ray eclipse # x (|Bradt fc McClintockl 
1 19831 ): 



1- I 0.38 -0.2 log — 



1/2 



(9) 



For systems with no measured orbital period, we make 
the standard LMXB assumption that md/m x < 0.8. Time- 
averaged accretion rate (as estimated from the X-ray flux) 
can give general constraints if we assume that mass transfer 
is driven by gravitational radiation: however magnetic brak- 
ing may also play a role in mass transfer and the contribution 
is hard to quantify. The conditions required for the system to 



iRappaport &: Jos 1 \ 19841) showed that helium-rich donor stars 
would have higher masses for a given P orD than expected for 
m ain sequence stars. However, the binary evolution simulations 
of lPodsiadlowski. Rappaport fc"P fahl ( 200^) did not generate any 
such stars at intermediate orbits, so we neglect this possibility un- 
less there is overwhelming evidence for the presence of helium-rich 
material (from X-ray bursts, for example). 



be persistent or t ransient at the inferred accre t ion rate were 
also considered dKing, Kolb fc Burderil Il996l ; iDubus et al.l 
1 19991 : [iri^t Zand. Jonker fc Markwardtll2007l ): unfortunately 
for most of the systems in our study this added very little 
in the way of tighter constraints. For systems where there 
are no better constraints on orbital period, we will assume 
that P or b lies between 10 minutes and 240 hours. The pre- 
sumed lack of eclipses in such systems sets i < 60° , so from 
Equation (JS| we obtain an upper limit on a x sini of 22.6 lt-s. 



3.3 Accreting millisecond pulsars 

The ten accreting millisecond pulsars are the only systems 
where we clearly know the NS spin frequency u s . All are 
transient, but their pulsation characteristics differ. Seven 
of these systems, which have short outbursts (weeks) and 
long periods of quiescence (years), show persistent pulsa- 
tions throughout their outburst phases. The other three sys- 
tems are rather different. HETE J1900. 1-2455 went into out- 
burst in 2005 and is still active, but only showed pulsations 
during the first two months of the outburst. The other two 
systems have shown pulsations only intermittently during 
outburst. We now discuss the sources in order of spin fre- 
quency, from highest to lowest. 



3.3.1 IGR J 0029 1+5934 (v 3 = 599 Hz) 

This source was discovered in 2004 during a 13 day out- 
burst with total fluenc e (0.1-200 keV) 1.8 x 10~ 3 ergs cm" 2 
(|Gallowav et all l2005t ). We assume that this is a typical 
outburst. The RXTE ASM history suggests a recurrence 
time of 3 years, giving P av = 1.8 x 10" 11 ergs cm" 2 s -1 . 
The most accurate position for the source, given by the op- 
tical/near infra-red counterpart, is R A = 00 h 29 m 03°05 ± 
(KOI, Dec = +59°34'18'.'93 ± C/05 (|Torres et alj 120081 ). 
This position accords with the earlier optical position of 
llFox fc Kulkarnill2004l ) and the Chandra X-ray (|Paizis et al.l 
120051 ) but is offset by 3.2<7 in RA from t he radio posi- 
tion |Rupen. Dhawan fc Mioduszewskil |2004| ). Mass trans- 
fer arguments sugges t a minimum distance of « 4 kpc 
(|Gallowav et alj 1200 5l ). and the lack of bur sts implies a 
maximum distance of 6 kpc (|Gallowavll2006l ). There have 
been three studies o f the spin and orbital parameters. 
iGallowav et"aH (|2005l ) used data from the start of the out- 
burst and the optic a l position of Fox fc Kulkarnil (|2004h . 
iFalanga et all (|2005l ): iBurderi et alj (|2007l ) used additional 
data from later in the outburst and the radio posit iorj 1 " 2 !. The 
studies agree on the orbital parameters: P or b = 8844.092 ± 
0.006 s, T, M = 53345.1875164 ± 4 x 10 -7 MJD (TDB), 
a x sini = 64.993 ± 0.002 It-ms (la uncertainties, Galloway, 
private communication), and e < 2 x 10" 4 (3a upper limit). 
There are however small differences in spin parameters. In 
our initi al assessment o f detec tability we use the most recent 
values of lBurderi et"aH (|2007l ): v s = 598.89213053±2 x 10" 8 
Hz, v a = (8.5±l.l)xl0 _13 Hz/sat epoch MJD 53346.184635 
(lcr uncertainties). 



12 Although an offset position can generate an apparent spin 
derivative, the offset is not sufficiently large to account for the 
reported values of v s - 
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3.3.2 Aql X-l (1908+005) (v s = 550 Hz) 

Aql X-l is a transient atoll source with quasi-regular out- 
b ursts, ~ 10 in the RXT E era. Using the spectral modelling 
of iGallowav et al.l (120081 1 we set 5 ASM cts/s = 3.1 x 10~ 9 
ergs cm -2 s _1 (bolometric). The long-term average ASM 
countrate is 1.8 cts/s, while the average countrate during the 
bright 45 day outburst of late 2000 was 34 cts/s. The most 
recent position, from VLA observations of the radio counter- 
part, is RA = 19 h li m lg;0 1±a;03, Dec = +00°35'06.'7±0'/3 
(J200 . la error bars ) jRupen, Mioduszewski fc Dhawanl 
120041 ). IGallowav et all (|2008l ) use PRE bursts to estimate 
a distance of 3. 9 ± 0.7 kpc o r 5.0 ± 0.9 kpc depending 
on composition. ICasella et all (|2008l ) have reported the 
detection of intermittent accretion-powered pulsations at 
550.2745 ± 0.0009 Hz. The orbital period, determined from 
the opt ical lightcurve in outburst, is P orb = 18.9479±0.0002 
hours (I Chevalier fc Ilovaiskvl 1 19981 ) (the quiescent period 
agrees at the 0.02% lev el). The tim e of m inimum opti- 
cal light determined by iGarcia et al.l (|l999l ) gives Tgo = 
2450282.220 ± 0.003 HJD. Both time of minimum light 
and orbital period were c onfirm e d in later analysis by 
IWelsh. Robinson fc Yound" (|200Ch . ICornelisse et al.l (|2007n 
attempted to measure the projected orbital velocity of 
the neutron star directly using phase-resolved optical spec- 
troscopy: they report a preliminary value of t> x sin i = 68 ± 5 
km/s. Although the fit quality is not good and this result re- 
quires confirmation, we will use this value in our preliminary 
assessment of detectability. 



3.3.3 SAX 1 1748.9-2021 (v B = 442 Hz) 

This is a transient source located in the globular cluster 
NGC 6440, with three o utbursts in the RXTE era. Using 
the spectral modelling of IGallowav et ah I (|200SD we set 13 
ASM cts/s = 5 x 10~ 9 ergs cm" 2 s" 1 (bolometric). The 
long-term average ASM countrate is 0.2 cts/s, while the 
average countrate during the bright 60 day outburst of 
2005 was 10 cts/s. T he position of the optical counterpart 
|Verbunt et al.l 200d) (which accords with the Chandra po- 
sition of lPoolev et al.l (|2002h l is RA = 17 h 4 8 m 52 ! :i4, Dec = 
-20°21 '32'/6 (J2000), with an error of 0'.'5. iKuulkers et~aH 
(|2003l ) estimate distance to the cluster as 8.4^J'g kpc, in 
good agreem ent with the value of 8 .1 ± 1.3 kpc derived from 
PRE bursts (|Gallowav et al.ll200St ). lAltamirano et al.1 (|2008l ) 
discovered intermittent accretion-po wered pulsations in tw o 
outbursts from this source (see also iGavriil et aD l|2007h ). 
IPatruno et al.l (|2008l ) have now carried out a detailed phase- 
connected timing study to determine the spin and orbital pa- 
rameters. They find v a = 442.36108118±5 x 10" 8 Hz, P orb = 



8.76525 ± 3 x 10~ J hours, T asc = 52191.507190 ± 4 x 10~ D 
MJD/TDB, and a x sini = 0.38760 ± 4 x 10~ 5 lt-s, with 
e < 1.3 x 10~ 4 (la uncertainties and upper limits). 

3.3.4 XTE J1751-305 (v a = 435 Hz) 

The 2004 outburst of this source, which lasted w 10 days, 
had total fluence (2-200 keV) (2.5 ± 0.5) x 10~ 3 ergs cm -2 
|Markwardt et al]|2002l ) . Assuming that this outburst is typ- 
ical, the mean recurrence time of 3.8 years yields P av = 
2 x 10~ n ergs cm~ 2 s~ 4 . The most accurate position for 
the source, measured with Chandra, is RA = 17 h 51 m 13!49± 



0!05, Dec = -30°37'23'.'4±0'.'6 (J20 00), where the uncertain - 
ties are the 90% confidence limits (jMarkwardt et al1l2002l ). 
Mass tra nsfer arguments yield a l ower limit on the distance 
of 6 kpc jMarkwardt et alll2002j ). The most recent timing 
study bv lPapitto et al. |200Sl ), using the Chandra position, 
gives P orb = 2545.342 ± 0.002 s, T asc = 52368.0129023 ± 
4 x 10~ 7 MJD (TDB), a x sini = 10.125 ± 0.005 lt-ms, e < 
1.3 x 10~ 3 (90% confidence level uncertainties and upper lim- 
its). The spin parameters are v a = 435.31799357 ± 4 x 10~ 8 
Hz, and the study also suggests that u a may be non-zero. 



3.3.5 SAX J1808. 4-3658 (v e = 401 Hz) 

IGallowav fc Cummind (|2006h analyse the five known out- 
bursts from this source, and find a mean outburst fluence 
of (6.0 ± 0.2) x 10" 3 ergs cm" 2 (0.1-200 keV). In estimat- 
ing transient gravitational wave signal we use an outburst 
duration of 20 days: we neglect the extended flaring phase 
often seen at the end of outbursts in this source, since emis- 
sion during this phase is at a much lower level. Mean re- 
currence time is 2.2 years, giving P av = 8.6 x 10 -11 ergs 
cm -2 s _1 . The most recent and precise position, for the op- 
tical counterpart, is RA = 18 h 08 m 27 i : 62, Dec = ^"SS^.^ 
(J2000), with an uncer t ainty of 0'.'15 (|Hartman et al.ll2008l ). 
IGallowav fc Cummind (|2006l ) have derived a distance of 3.4- 
3.6 kpc using both mass transfer arguments and burst prop- 
erties. Using the refined optical position, lHartman et al.l 
(|2008l ) report the following values for the orbital param- 
eters: Porb = 7249.156961 ± 1.4 x 10" 5 s at time T asc = 
52499.9602477 ± 1.0 x 10~ 6 MJD (TDB), P orb = (3.48 ± 
0.12) x 10~ 12 Hz/s, a x sini = 60.28132 ± 2.4 x 10~ 4 lt-ms 
(la errors). The eccentricity e < 0.00021 (95% upper limit). 
The spin rate was tracked across multiple outbursts, and is 
given relative to a reference frequency vq — 400.975210 Hz. 
In the 1998 outburst, v s - i/ = 0.371 ± 0.018 ^Hz, with v s 
in the range (-7.5, 7.3) xl0~ 14 Hz/s (95% confidence lim- 
its). In the 2000 outburst, v s - u = 0.254 ± 0.012 ^Hz, 
with v a in the range (-1.1, 4.2) xl0~ 14 Hz/s. In the 2002 
outburst, v a — i/o = 0.221 ± 0.006 /iHz, with i> a in the range 
(-1.2, 2.5) xl0~ 14 Hz/s. In the 2005 outburst, v s - v = 
0.190 ± 0.015 ^Hz with v a in the range (-0.5, 2.3) xlO" 14 
Hz/s. Fitting the frequency evolution across all four out- 
bursts gives i> a = (5.6 ± 2.0) x 10~ 15 Hz/s. In our initial 
assessment of detectability, we use the spin solution for the 
2002 outburst. 



3.3.6 HETE .11900.1-2455 (v a = 377 Hz) 

This source was first detected in June 2005, and ha s 
remained in ou tbur st ever since (iDegenaar et all 120071 ). 
iGallowavl (|2006l ) and IGallowav et al. I (|2007al ) report an av 
erage (bolometric) outburst flux of ~ 9 x 10 



10 ergs cm 2 



s . We assume an outburst duration of 2 years, and (given 
that no previous out bursts are known) a recurrence time 
of at least 10 years (|Kaaret et al.l 120061 ). yielding P av w 
2 x 10~ 10 ergs cm -2 s~ 4 . The position of the optical coun- 
terpart, is RA = 19 h 00 m 08"65, Dec = -24° 55' 13? 7 (J2000), 
with an estimated uncertainty of d .'2 (Fox 2005). The dis- 
tance estimated by IGallowav et ah I (|2008l) using RXTE ob- 
servations of radius expansion bursts, is 4.7 ± 0.6 kpc. 
This accords with the earlier estimate of 5 kpc made by 
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iKawai fc Suzukil (l2005h using H ETE burst data. Timing 
analysis by Kaaret et alj (|2006l ). using RXTE data from 
June 16 - July 7 2006, resulted in the following orbital pa- 
rameters: P orb = 4995.258 ± 0.005 s, T 90 = 53549.145385 ± 
7x 10~ 6 MJD (TT), a x sini = 18.41 ±0.01 lt-ms, e < 0.002. 
The spin was u s = 377.296171971 ± 5 x 10~ 9 Hz. All errors 
and upper limits are la uncertainties. On July 8th there 
was an apparent jump in spin rate of Av s /v s ~ 6 x 10~ 7 , 
to v s = 377.291596 ± 1.6 x 10 -5 Hz. Thereafter pulsations 
ceased and spin has not been tracked since despite the fact 
that t he source has remained in outburst ('Ga lloway et al.l 
l2007al ). In our initial assessment of detectabihty we do not 
take into account the apparent jump in spin. 

3.3.7 XTE J1814-338 (v„ = 314 Hz) 

This pulsar has ha d only on e outb urst in the RXTE era, 
lasting « 50 days. IGallowavl (|2006l ) estimate a bolometric 
outburst fluence of (3.0 ±0.1) x 10 -3 ergs cm -2 s _1 . Given 
a recurrence time of at least 7.5 years, this yields an up- 
per limit on P av of 1.3 x 10 -11 ergs cm -2 s _1 . The posi- 
tion derived from X-ray and optical spectroscopy is RA = 
18 h 13 m 39:04, Dec = -33°46'22"3 (J200 0), 90 % confidence 
error circle of 0'.'2 (|Krauss et alj [20 051). These authors use 
X-ray spectroscopy to infer a minimum distance of 3.8 kpc. 
The upper limit on the distance, derived fro m burst proper- 
ties, i s 8.0±1.6 kpc (Strohma yer et ai1 l2003f). Timing analy- 
sis bv lPapitto et al. I (|2007h . using data from the whole 2003 
outburst, leads to the following orbital parameters: P or b = 
15388.7229±0.0002 s, T asc = 52797.8101698±9 x 10" 7 MJD 
(TDB), a x sini = 390.633 ± 0.009 lt-ms, e < 2.4 x 10~ 5 
(3a upper limit). The associated spin parameters are v 3 = 
314.35610879± 1 x 10~ 8 Hz, v s = (-6.7±0.7) x 10~ 14 Hz/s. 
Uncertainties are at the 90% confidence level. 



3.3.8 XTE J1807-294 (v B = 191 Hz) 

There ha s been only one r ecorded outburst from this source, 
for which IGallowavl (|2006h computes a bolometric fluence of 
(3.1 ± 0.2) x 10~ 3 ergs cm -2 . For a recurrence time of at 
least 7.1 years, F av is at most 1.4 x 10 -11 ergs cm -2 s _1 . 
In estimating transient gravitational wave signal we use an 
outburst duration of 50 days, neglecting the prolonged low 
flux tail at the end of the outburst. The most accurate posi- 
tion, measured with Chandra, is RA = 18 h 06 m 59:8, Dec = 
-29° 24' 30" ( J2000), with an uncertainty due to systematic 
errors of 1" dMarkwardt, Juda fc Swankll2003|) . Using mass 
transfer estimates. IGallowavl ( 2006h derives a lower limit on 
the distance of 4.7 kpc. Determination of spin and orbital 
parameters in this source is complicate d by extreme varia- 
tions in pulse profile. A recent study bv lRiggio et al.1 (|2007l ) 
finds P orb = 2404.41665 ± 0.00040s, a x sin i = 4.819 ± 0.004 



3.3.9 XTE J0929-314 (u s = 185 Hz) 



lt-ms, T asc = 52720.675603 ± 6 x 10~ b MJD/TDB and 
v s = 190.62350694 ± 8 x 10 8 Hz (la uncertainties ) The 
eccentricity e < 0.0036 (2a upper limit). IChou et al.1 (|2008l ) 
report similar values apart from for r asc , where the values 
found by the tw o studies differ by m ore than the quoted un - 
certainties. Both lChou et al.1 (|2008l ) and lRiggio et al.1 (|2008l ) 



suggest a non-zero i> s 



10" 



10" 



Hz/s in outburst. 



This s ource has had o ne outburst in the RXTE era, for 
which IGallowavl (|2006l ) estimates a bolometric fluence of 
(5.4 ± 0.3) x 10 -3 ergs cm -2 . Outburst duration, which we 
assume to be typical, was ~ 60 days. Given a recurrence time 
of at least 6.3 years, F av is at most 2.7 x 10 -11 ergs cm" 2 
s - . The most accurate position, given by the optical coun- 
terpart, is RA = 9 h 29 m 20:i9, Dec = - 31 o 23'03'.'2 ( J2000) 
with error circle O'.'l |Giles et alj|2005h . IGallowavl (|2006l ) 
uses mass transfer and recurrence time estimates to infer 
a lower limit to the distance of 3.6 kp c (revising an earlier 
value of 5 kpc in Galloway et al.l ((2002)). Timing analysis by 



In our initial assess ment of detectabihty we use the values 
and uncertainties of iRiggio et all (|2007n . although clearly 
the 'true' uncertainty on T asc is larger. 



Galloway et al.l (120021) . using an earlier optical position from 
Giles et al.1 (|2002l ~ leads to the following orbital parameters: 
P orb = 2614.746 ± 0.003 s, T 90 = 52405.49434 ± 1 x 10 -5 
MJD (TDB), a x sini = 6.290 ± 0.009 lt-ms. The associated 
spin parameters are v s = 185.105254297 ± 9 x 10~ 9 Hz, 
v s = (-9.2 ± 0.4) x 10~ 14 Hz/s. All errors are la uncer- 
tainties (Galloway, private communication). The eccentric- 
ity e < 0.007 (2a limit). 



3.3.10 SWIFT J 1756. 9-2508 (v 3 = 182 Hz) 

This source has had one outburst in 2006 lasting 13 days, 
with total fluence (1-10000 keV) (4.5±0.8) x 10~ 4 ergs cm" 2 . 
No previous outbursts are known (although there are gaps 
in coverage), but this suggests a recurrence time of at least 
10 years. The best position for the source, from SWIFT, is 
RA = 17 h 56 m 57:35, Dec = -25°06'27'.'8 (J2000), with un- 
certainty 3'.'5. The distance is not well constrained, but is 
thought to be ~ 8 kpc. The orbital parameters are P or b = 
3282.104 ± 0.083 s, T asc = 54265.28707 ± 6 x 10~ 5 MJD 
(TDB), a x sini = 5.942 ±0.027 lt-ms, e < 0.026 (95% upper 
limit). The spin rate is v s = 182.065804253 ± 7.2 x 10~ 8 Hz, 
with v s < 1 x 10~ 12 Hz/s. All errors are 90% confidence, up- 
per limits are 95% confidence. All information on this source 
is taken from lKrimm et all (|2007l ). 



3.4 Burst oscillation sources 

We list all sources for which burst oscillations have been 
reported, in order of decreasing burst oscillation frequency 
vi, . For some sources, oscillations have been detected at the 
same frequency in multiple bursts: these results can be re- 
garded as secure. In some cases, however, oscillations have 
only been seen in a single burst. Given the number of NSs 
whose bursts have now been searched for oscillations (a fac- 
tor not included in quoted statistical significances), these 
results should be regarded as tentative until confirmed in a 
second burst. We include them in our survey (marked with 
an asterisk) since if they turn out to be promising gravi- 
tational wave sources this would provide added impetus to 
confirm or alternatively rule out the candidate burst oscil- 
lation detection. 



3.4.1 XTE J 1739- 285 (u b = 1122 Hz)* 

This transient atoll source has had four out bursts in the 
RXTE era. Using the spectral modelling of iKaaret et al.l 
|2007l ) we set 10 ASM cts/s = 1.4 x 10~ 9 (2-20 keV) and ap- 
ply a bolometric correction factor of 1.34 (the average factor 
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found bv lGallowav et aT] (|2008l ) for burst sources) . The long- 
term average ASM countrate is 0.3 cts/s, and the average 
countrate during the bright 200 day outburst of 2005 was 4 
cts/s. The most precise source position, measured by Chan- 
dra, is RA = 17 h 39 m 5ff!95, Dec = -28°29'46"8 (J2000), with 
a 90% error radius of 0'.'6 ijKrauss et alj|200d ). The absence 
of an optic al/IR counterpart s ets an upper limit on distance 
of 12 kpc (|Torres et al.l 120061 '). The absence o f PRE bursts 
sets a more stringent upper limit of 10.6 kpc l)Kaaret et al.l 
l2007h . A candidate 1122 Hz burst oscillation was detected in 
part of one burst recorded by RXTE during the 2005 out- 
burst. The orbital parameters of the source are unknown, 
and there are no measured constraints on the properties of 
the companion. 

3.4.2 4U 1608-522 (Ub = 620 Hz) 

This transient atoll source has had several outbursts dur- 
ing RXTE's lifetime, and seems to be active at a low level 
(ASM countrate > 3a) e ven when not in outb urst. Using 
the spectral modelling of iGallowav et ah I (|2008h we set 70 
ASM cts/s = 4.8 x 10~ 8 ergs cm" 2 s" 1 (bolometric). The 
long-term average ASM countrate is 3.6 cts/s, while the av- 
erage countrate during the bright 100 day outburst in 2005 
was 30 cts/s. The best position, from the optical counter- 
part, is RA = 16 h 08 m 52 i ;2, Dec = -52° 17'43" (B1950), with 
errors of ±0'.'5 (|Grindlav fc LiUedll978T ). No more up to date 
position is available. Source distance, derived under the as- 
sumption that radius expansion bursts rea ch the Edding- 
ton l imit for pure helium, is 4.1 ± 0.4 kpc (|Gallowav et al.l 
I2008T I. Burst os cillations at ~ 620 Hz have been detected i n 
multiple bursts (|Hartman et al . 20031 IGallowav et alJliooj ). 
The last (unsuccessful) search for persistent pulsations from 
this source, using data from 198 9-1991, searched only up 
to 512 Hz jVaughan et al.lll994l ). The orbital parameters 
are not known, although a number of tentative periodic- 
ities have been rep orted (jLochner fc Roussel-Duprel 1 19941 ; 
IWachter et ai]|2002l ). The spectral type of the companion is 
that of a late F/early G type main sequence companion star, 
but would also match that of a more evolved K/M type star 
jWachter et al.ll2002l ). The mean density range for donors of 
this type suggests that the orbital period lies in the range 
10 - 125 hours. Given that the system is non-eclipsing we 
assume i < 60° and hence obtain, for the assumed orbital 
periods, a x sini < 12.3 lt-s. 

3.4.3 SAX J1750. 8-2900 (v h = 601 Hz) 

This weak transient atoll source has had two outbursts 
in the RXTE era. F ollowing the spectral modelling of 
IGallowav et~aT] (|2008f ) we set 8 ASM cts/s = 3.4 x 10~ 9 
ergs cm" 2 s" 1 (bolometric). The long-term average ASM 
countrate is 0.08 cts/s, while the average countrate during 
the brighter 110 day outburst in 1997 was 2.8 cts/s. The 
most precise position, from BeppoSAX, is RA = 17 h 50 m 24 s , 
Dec = -29° 02' 18" (J2 000), with a 99% error radius of 1' 
l|Natalucci et al.ll 19991 ). The distance, estimated from radiu s 
expansion bursts, is 6.79 ± 0.14 kpc (|Gallowav et alj|2008h . 
Burst oscill ations at 601 Hz have been detected in multi- 
ple bursts (|Kaaret et alj|2002l : IGallowav et all 120081 ). The 
orbital parameters are not known, and there are no obser- 
vations of the companion star. 



3.44 GRS 1741.9-2853 (AX J1745. 0-2855) (u b = 589 
Hz) 

This is a transient source in the crowded Galactic Centre, 
with three outbursts in the RXTE era. Unfortunately the 
fact that the field is crowded means that there is no reliable 
ASM flux history for this source. The most precise position 
for the source, derived by Chandra, is RA = 17h45 m 2 i !33, 
Dec = -28°54'49"7 (J2000) with an uncertainty of 0'.'7 
l|Muno. Baganoff fc Araba diis 2003|). Analysis of radius ex- 
pansion bursts suggests a distance 6.0 ± 1.6 kpc (Eddington 
limit assuming cos mic abundances) or 8 ± 2 kpc (Eddington 
limit for pure He) (|Gallowav et al.l 120081 ). Burst oscillations 
at 58 9 Hz have been detected in 2 bursts |Strohmaver et aT] 
19970 The orbital parameters are not known, and there is 
no information on the properties of the companion star. 



3.4.5 4U 1636-536 (v b = 581 Hz) 

This is a persistent atol l sourc e. Using the spectral mod- 
elling of IGallowav et~ai1 (|2008l ) we set 18 ASM cts/s = 
8.4 x 10 -9 ergs cm -2 s _1 (bolometric). The long-term av- 
erage ASM countrate is 10 cts/s. The best position for 
the source, from the optical counterpart V801 Ara, is RA 
= 16 h 36 m 56":41, Dec = -5 3°39'18'/1 (B1950), with a n er- 
ror circle of less than 1" (|Bradt fc McClintock1ll983l) . No 
more up to date position is available. The distance, esti- 
mated from a large sample of radiu s expansion bursts, is 
6.0 ± 0.5 kpc (jGallowav et al.1 12006? ). Burst oscillations at 
« 581 Hz have been seen in multiple bursts and a super- 
burst (IStrohmaver et all Il998al ; IStrohmaver fc Markwardtl 
120021 ; IGallowav et al.l 20081 ). The most recent unsuccess- 
ful search for persistent pulsat ions, using data from 1987, 
searched only up to 512 Hz (|Vaughan et al.1 Il994l ). The 
latest ephemeris, derived from phase-resolved optical spec- 
troscopy, gives P OI .b = 0.15804693 ± 1.6 x 10" 7 da ys an d 
T 90 = 24 528 1 3 531 ± 0.002 HJD (|Casares et all ' I2OO6I ). 
lAugusteiin et al.l (|l998l ) set an upper limit on the orbital 
period derivative of | -Porb /-Port. | ^ 3 x 10 s years. No tighter 
limit has yet been reported. Orbital Doppler shifts on burst 
oscillations during a superburst lead to limits 90 < « x sin i < 
113 k m/s, l|Strohmaver fc Markwardtl |2002| ; ICasares et al.l 
120061 ). The quoted range corresponds to varying the refer- 
ence phase of the ephemeris across the ±la range. 



3.4.6 X 1658-298 (MXB 1659-29) (v b = 567 Hz) 

This transient source, which has eclipses and dips, was ac- 
tive from 1976-1979, and then again from 199 9- 2001 Us- 
ing the spectral modelling of IGallowav et al.l (|2008l ) we 
set 2.4 ASM cts/s = 6.7 x 10" 9 ergs cm" 2 s" 1 (bolo- 
metric). Over the most recent outburst, which lasted 870 
days, the average ASM countrate was 2.7 cts/s. Assum- 
ing a recurrence time of « 23 years, the long-term aver- 
age ASM countrate is 0.3 cts/s. The best position for the 
source, from the optical counterpart V2134 Oph, is RA = 
17 h 02 m 06?42, Dec = -29°56 / 44 / /33 (J2 000), with an uncer- 
tainty of O'.'l (|Wachter fc Smald UJoj ) . Assuming that the 



13 The bursts were originally attributed to a different source, 
MXB 1743-29, due to source confusion in the crowded field 
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bright ra dius expansion bursts r each the Eddington limit for 
pure He, iGallowav et al.l (|2008l ) derive a distance of 12 ± 3 
kpc. Burst oscillat ions at ~ 567 Hz have been detecte d 
in multiple bursts |Wiinands. Strohmaver fc Franco! 200ll). 
The m ost up to date X-ray ephemeris, by Qosterbroek et all 
(|200ll ). gives P orb = 0.2965045746 ± 3.4 x 10" 9 days, 
and T 90 = 2443059.225826 ± 0.000093 JD/TDB. The pro- 
jected semi-major axis has not been measured, but we 
can set bounds on it using the constraints outlined at 
the start of this section . For the above orbital period, 
IWachter. Smale fc Baiiynl \200(t ) find X-ray eclipse dura- 
tion half-an gles 6.34° ± 0.01° (sl i ghtly lower than the value 
reported by 



les 6.34 ± U.U1 (sl i ghtly 1c 
Cominskv fc Wood! (|l984l )) 



. We can therefore 
we assume Roche 



use Equation to restrict sin i. If 
lobe overflow, then from Equation (Q the mean density 
of the donor star is 2.2 g/cm 3 , suggesting a main se- 
quence or evolved companion. Maximum donor mass oc- 
curs if the donor is on the main sequence. In this case, 
max imum Mg = 0.78Mp) a nd would be of spectral type 
K0 (|Wachter fc Smalei Il998h . Whilst the data are consis- 
tent with this spectral type, there is some indication that 
the spectral type is later, suggesting a more evolved (lower 
mass) companion with spectral type perhaps as late as M2 
|Wachter. Smale fc Bailvn]|2000l ). We will therefore adopt a 
minimum c ompanio n mass of 0.1M ? jBaraffe fc Kolbll2000l ; 



IPodsiadlowski. Rappaport fc Pfahlll2002l) . For the neutron 
star we will consider ma sses in the range 1.2 — 2.4M© 
IlLattimer fc PrakashlbOOTf ). These assumptions suggest that 
dxSini lies in the range 0.24-2.03 lt-s. 

3.4.7 A 1744-361 (v b = 530 Hz)* 

This transient dipping source has irregular outbursts, three 
in RXTE's lifetime. Based on the spectral modelling of 
iBhattacharvva et all (|2006l ) we set 1.2 ASM cts/s = 3 x 
10~ 10 ergs cm~ 2 s _1 (3 - 14 keV) flux and apply a bolo- 
metric correction factor of 1.34 (the mean correction factor 
for converting 2 . 5 - 25 keV flux to bolometric flux found by 
IGallowav et al.l (|2008l ) for other burst sources). The long- 
term average ASM countrate is 0.16 cts/s, and during the 
bright 100 day outburst in 2003 the average countrate was 
6 cts/s. The best position, for the radio counterpart, is RA 
= 17 h 48 m 13" 148 ±0^014, Dec = -36°07 , 57 / /02 ± 0"3 (J2000) 
ijRupen. Dhawan fc Mioduszew ski 2003| ). This is within th e 



error circle of the Chandra position (|Torres et al.l l2004f ) . 
The lack of radius expansion in the one burst detected 
by RXTE sets an upper li mit on the distance of « 9 kpc 
|Bhattacharvva et al.1120061 ). These authors reported a burst 
oscillation at ~ 530 Hz during the rising phase of the one 
burst detected by RXTE. They also reported a possible or- 
bital period of 97 ± 22 minutes, traced by dips in the X-ray 
lightcurve, (jBhattacharvva et al.1120061) . In our initial assess- 
ment we will assume that this is indeed the orbital period. In 
this case the assumption of Roche lobe overflow (Equation 
[7]) gives a mean donor density of 28 — 70 g/cm 3 , consistent 
with a main sequence star or slightly evolved donor. The 
maximum donor mass (for a star on the main sequence) is 
0.22 M®. For a dipping source we can assume an inclination 
in the range 60 — 75° . For neutron star masses in the range 
1.2 — 2AMq we can therefore set an upper limit on a x sini 
of 0.2 lt-s. The minimum donor mass, of 0.07Mq, is set by 
the most evolved hydrogen burning star possible. Together 



with the inclination constraint, this gives a lower limit on 
a x sin i of 0.05 lt-s. 



3.4.8 KS 1731-260 (u b = 524 Hz) 

This transient atoll source was in outburst from 1988 (and 
p ossibly earlier) until 20 01. Based on the spectral modelling 
of IGallowav et all (|2008h we set 18 ASM cts/s = 9.7 x 10~ 9 
ergs cm -2 s -1 (bolometric). The average countrate while 
the source was still in outburst was 9 cts/s: the long-term 
average countrate over the RXTE lifetime is 4 cts/s. The 
most precise position, for the optical counterpart, is RA 
= 17 h 34 m 13 ! !47, Dec = -26°05'18'.'8, with an error of 0'.'4 
|Wiinands et al.ll2q0l|). This is within the error circle of the 
Chandra position (|Revnivtsev fc Sunvaevll2002l ). Assuming 
tha t the radius expansio n bursts reach the limit for pure 
He. IGallowav et al.l (|2008h derive a distance of 7.2 ± 1.0 kpc. 
Burst oscil lations at ~ 524 Hz have been detected in multi- 
ple bursts (ISmith. Morgan fc Bradlll997l ; iMuno et al1l2000l ; 
IGallowav et al.ll200Sft . The orbital parameters have not been 
measured. A study o f the scatter of asy mptotic burst oscil- 
lation frequencies bv lMuno et al.l (|2000l ) suggested that the 
v x sin i might be as high as 340 ± 100 km/s, but this has not 
been revisited. Identification of the counterpart was ham- 
pered by high redde ning along the galact ic plane, but it has 
now been detected jMignani et al.ll2002l ). If the companion 
is on the main sequence, it has to be of spectral type later 
than F: if it has evolved off the main sequence then it is not a 
red giant. For this to be the case in a Roche lobe overflowing 
system we require P or b > 2 hours. 



3.4.9 4U 0614+09 (v b = 415 Hz)* 

This is a bursting atoll source, persiste nt but high l y var i- 
able. Based on spectral modelling by I Ford et "all (|2000t ). 
we set 9 ASM cts/s = 3.3 x 10~ 9 ergs cm -2 s" 1 
(bolometric). The long-term average ASM countrate is 
3.3 cts/s. The best and most recent position, from a 
Spitzer observation of the IR counterpart, is RA = 
6 h 17 m 07:35 ± 0103 Dec = +09°08'1 3 .'60 ± 0'.'05 (J2000) 
|Migliari et al.l |2006| PI iBrandt et all l|l992f ) infer an up- 
per limit on the distance of 3 kpc from an X-ray burst. 
IStrohmaver. Markwardt fc Kuulkersl |2008l) detected burst 
oscillations at 415 Hz in one burst recorded by the 
SWIFT Burst Alert Telescope. The orbital parameters 
are u nknown, but this is a candidate ul tra-compact bi- 
nary (|juett. Psaltis fc Chakrabarty||200fl ). iNelemans et al.l 
|2004l ) have shown that the companion is most likely 

a C/O white dw arf. Using the white dwarf models of 

iDelove fc Bildstenl (|2003h this would imply P orb = 15-20 
minutes and hence (assuming that i < 60° due to the lack 
of dips and eclipses), a x sini ^ 0.014 lt-s. 



3.4.10 4U 1728-34 (GX 354 +00) (u b = 363 Hz) 

This is a p ersistent atoll s o urce. Based on the spectral 
modelling of IGallowav et al.l (|2008l ) we set 4 ASM cts/s = 



14 This is consiste nt with the older position fo r the optical coun- 
terpart V1055 Ori llBradt fc McClintocklll983h . 
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1.2 x 10 -9 ergs cm~ 2 s _1 (bolometric) . The long-term av- 
erage ASM countrate is 7.3 cts/s. The most precise posi- 
tion for this source is that of the radio counterpart, RA 
= 17 h 31 m 57:73 ± 0!02, Dec = -33°50 / 2 / /5 ± l'.'l (J2000) , 
where the errors are la uncertainties (|Martf et ai.l 1 199a ) ■ 
Assuming that the b right PRE bursts reach the He limit, 
iGallowav et alj (|2008l ) infer a distance of 5.2±0.5 kpc. Burst 
osci llations have been detec t ed in multiple bursts at « 363 
Hz (jStrohmaver et al.lll996l : IGallowav et al.ll200St ). An un- 
successful search fo r persistent pulsations using Ginga data 
was carried out by IVaughan et al (Il994 ). The orbital pa - 
rameters are not known. although lStrohmaver et al 1 (|l998bh 
infer u x sin i < 20.7 km/s from the scatter of asymptotic fre- 
quencies of burst oscillations from a series of bursts from 
1996-7. We will use this constraint in our initial assessment 
and use it to estimate P OT b- If we assume Roche lobe over- 
flow, then we find that a main sequence or evolved star can- 
not satisfy the various relations. The companion must be ei- 
ther a white dwarf or a helium star. This is consistent with 
the properties of the X-ray burs ts from this source, w hich 
suggest a hydrogen poor donor (jGallowav et al.|[2008l ). We 
will therefore assume that P orb < 10 hours. 



34.ll 4U 1702-429 (v b = 329 Hz) 

This is a p ersistent atoll source. Based on the spectral 
modelling of IGallowav et all l|2008t) . we set 2 ASM cts/s = 
7.8 x 10 -10 ergs cm~ 2 s (bolometric). The long-term aver- 
age ASM countrate is 3.2 cts/s. The best position, measured 
by Chandra, is RA = 17 h 06 m 15!314 Dec = -43°02'08'.'6 
(J2000), with an uncertainty of 0'.'6 ijWachter et al.ll2005h . 
Assuming that t he PR E bursts reach the limit for pure He, 
IGallowav et alj (|2008l ) derive a distance of 5.5 ± 0.2 kpc. 
Burst oscillations are detected in multiple b u rsts at 329 Hz 
llMarkwardt. Strohmaver fc Swankl Il999al : IGallowav et al.l 
The orbital parameters are not known but there are 
no dips or eclipses. 



3.4.12 MXB 1730-335 (Rapid Burster) (v b = 306 Hz)* 

The Rapid Burster is a transient globular cluster source with 
regular outbursts that occur around every 200 days. It is un- 
usual in being the only system to show both Type I and Type 
II X-ray bursts, the latter being driv en by spasmodic accre- 
tion. Using the spectral modelling of IGallowav et all (|2008l ) 
we set 10 ASM cts/s = 6.7 x 10-9 ergs cm -2 s" 1 (bolo- 
metric). The long-term average ASM countrate is 1 ct/s, 
and during a typical 25 day outburst the average countrate 
is 8.8 cts/s. The most accurate position, given by the ra- 
dio counterpart, is RA = 17 h 3 3 m 24°61, Dec . = -33 °23 / 19"8 
(J2000), with an error of O'.'l (|Moore et alj|2000h. This is 
withi n the error circle of Chandra observations ( Homer et al.l 
l200ll ). The distance to the host g lobular cluster, Liller 1, is 
S.Siji kpc (|Kuulkers et alj|2003 ). Averaging the burst rise 
phase of 31 X-ray bursts recorded by RXTE revea led a weak 
cand idate burst oscillation frequency of « 306 Hz (|Fox et al.1 
l200lf ). The orbital parameters are not known, and no opti- 
cal counterpart has been detected because of crowding in 
the host globular cluster. 



3.4.13 IGR J17191-2821 (u b = 294 Hz)* 

This transient X-ray source was discovered only recently, 
with one recorded outburst. The average ASM countrate 
during the 11 day outburst was 3.7 cts/s, giving a long-term 
av erage countrate du r ing the RXTE era of 0.01 cts/s. Based 
on iKlein-Wolt et al.1 (|2007aD we assume 3.5 ASM cts/s = 
1.2 x 10" 9 ergs cm" 2 s" 1 (2-10 keV), and apply a bolo- 
metric correction factor of 2. The most precise position for 
the source, from SWIFT, is RA = 259. 81306°, Dec = - 
28.299 19° (J2000), with an accuracy of 4" (|Klein-Wolt et al.1 
l2007bh . Burst oscillations at fti 294 Hz have been detected 
in one burst (|Markwardt et al.1 [20071 ). The peak flux of the 
X-ra y bursts sets an upper limit to the distance of ~ 11 
kpc (jMarkwardt et alll2007l ) ■ Orbital parameters are not yet 
known, and there is no information on the companion. 



3.4.14 4U 1916-053 (u b = 270 Hz)* 

This is a persistent source in an ultracompact binary 
with an H-poor donor s t ar. F ollowing the spectral mod- 
elling of IGallowav et""ai1 (120081 ) we set 0.5 ASM cts/s = 

2.7 x 10 -10 ergs cm -2 s _1 (bolometric). The long-term av- 
erage ASM countrate is 1.3 cts/s. The most recent posi- 
tion, given by Chandra, is RA = 19 h 18 m 47:871, Dec = - 
05°14 / 17'.'09 (J2000) with an error 0'.'6 (|laria et allbOOrt Fl. 
IGallowav et al.1 (|2008l ) derive a distance of 8.9 ± 1.3 kpc or 

6.8 ± 1.0 kpc, depending on composition, from PRE bursts. 
Burst oscilla tions at m 270 Hz ha ve been detected in one 
X-ray burst (|Gallowav et al.l|200l| ). The orbital period has 
been the subject of much debate in the literature, due in 
part to differences between the X-ray and orbital periods 
determined by dipping and photometry respectively. The 
most recent papers on this topic seem to resolve t he issue 
jChou. Grindlav fc Bloserl l200ll ; iRetter et all |2002| ) by de- 
termining that the system displays superhumps rather than 
being a hierarchical triple. IChou. Grindlav fc Bloserl (|200ll ) 
use X-ray dip times to derive P or b = 3000.6508 ± 0.0009 s 
(in agreement with the v a lue of P olh = 3000.6452 ± 0.0043 
s reported by IWen et all (|2006l ) using the RXTE All Sky 
Monitor), with P orb < 2.06 x 10" 11 (2a upper limit). The 
X-ray dip ephemeris is T di p = 50123.00944±1.4x 10" 4 MJD. 
Whilst it is not entirely clear how the dip time relates to T90 
we will assume in our initial analysis that this could be de- 
termined, and use the uncertainty in Tdi p as an estimate of 
the uncertainty in T90. The projected semi-major axis has 
not been measured, but we can impose some constraints. 
For the known orbital period, the assumption of Roche lobe 
overflow (Equation [7} gives the mean mass of the donor 
as 158 g/cm 3 , implying a dwarf companion. X-ray burst 
properties imply a helium rich donor. Optical spectroscopy 
shows large amounts of N, suggesting that the companion 
is a helium white dwarf rather than a helium star or an 
evolved secondary (iNelemans. Jonker fc Steeghsll2006l ). Us- 
ing the models of iDelove fc Bildstenl (|20ol T this implies a 
companion mass in the range 0.008 - 0.03 M© (depending 
on core temperature). For dips we expect an inclination in 
the range 60 — 75° . If we assume that the neutron star mass 



15 This position has a substantial offset in decli nation from the 
optic al counterpart (V1405 Aql) reported by llGottwald et all 
Il99ll) . with RA = 19 h 18 m 47!91, Dec = -05°14'08'.'7 (J2000). 
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is in the range 1 .2 — 2AMq , Equation (J5J implies that a x sin i 
must lie in the range 4-25 lt-ms. 

3.4.15 XB 1254-690 (u b = 95 Hz)* 

This source is persi stent, at a steady level, with dips. 
I Galloway et all (|2008l ) carry out spectral modelling and find 
a persistent flux level of 9 x 10~ 10 erg c m" 2 s -1 (bolomet - 
ric). The Chandra position reported by llaria et alj (|2007l ) 
is RA = 12 h 57 m 37tl53, Dec = -69° 17' 18'.' 98, with a 90% 
uncertainty radius of tf/rF^l Radius expansion in the precur- 
sor t o a superburst leads t o a distance est i mate of 13 ± 3 
kpc (jin't Zand et alj|2003h . iBhattacharvval (|2007i ) reported 
tentative evidence of a burst oscil l ation at 95 Hz in one 
burst from this source. iMotch et all (|l987l ) used optical pho- 
tometry to find Porb = 3.9334 ± 0.0002 hours This accords 
with the period of 3.88 ± 0.15 hours d erived from the X- 
ray lightcurve ()Courvoisier et all 198rJ). although the X - 
ray dips are not al ways present (jSmale fc Wachterl Il999l ) . 
iBarnes et al. (I2007T) ha ve rec ently updated the ephemeris 
derived bv lMotch et al.l ij 19871 ') and report an X-ray dip time 
T di p = 2453151.647 ± 0.003. iMotch et all ij 19871 ) showed 
that the X-ray dips in this source occurred at phase 0.84 
(with zero at optical minim um), so they therefo re derive To 
= JD 2453151.509 ±0.003. IBarnes et all (|2007m have used 
phase-resolved spectroscopy of the He II A4686 emission line 
(thought to be emitted in the inner accretion disk, close to 
the compact object) to estimate the velocity of the compact 
object, and find a velocity semi-amplitude u x sini = 130±16 
km/s. Only the lower portion of this range is consistent with 
a main sequence or undermassive companion star: however 
we will use this value in our initial assessment of detectabil- 
ity. 



S.4.I6 EXO 0748-676 (v b = 45 Hz) 

This system, which shows both dips and eclipses, has been 
persisten tly active since 1 985. Based on the spectral mod- 
elling of iGallowav et all (|2008l ), we set 0.6 ASM cts/s = 
3.6 x 10~ 10 ergs cm -2 s _1 (bolometric). The long-term aver- 
age countrate is 0.76 cts/s. The best position, for the optical 
counterpart UY Volantis, is RA = 7 h 48" 1 25°0 ± 0?1, Dec = 
-67°37'31'.'7±0'.'7 (B1950) (|Wade et alJll985l ). No more up 
to date position is available. The detection of PRE bursts 
that seem to be He ri c h implies a distance of 7. 4 rfc 0.9 kpc 
(jGallowav et al.ll2008l ). IVillarreal fc Strohmaveil <|2004h dis- 
covered burst oscillations at 45 Hz by averaging together 
spectra from 38 separate bursts detected between 1996 and 
early 2003. Timing of the eclipses constrains the orbital pa- 
rameters. Attempts to compute an orbital ephemeris, how- 
ever, have been complicated. The most recent study by 
IWolff et alj (|2002l ) finds a large apparent period change of 
8 ms over the period 1985 to 2000 - much larger than ex- 
pected from orbital models - and intrinsic jitter that cannot 
be explained by any simple ephemeris. The reason for this 
variability has yet to be resolved, and it is not clear whether 
this represents genuine evolution in the binary period or not. 



16 This position is 2" away from the older opt ical counterpart 
position reported bv lBradt fc McClintockl jl983ft . 



IWolff et al l (|2002l ) consider various models in their analy- 
sis, with orbital periods P or b = 0.1593378 ± 1 x 10~ 7 days, 
mid-eclipse times T, M = 46111.0739 ± 0.0013 MJD/TDB 
and I Porb I < 10 -11 . More recent analysis by IWen et al] 
(|2006r i. using ASM data, finds P orb = 0.1593375 ± 6 x 10" 7 
days. The projected semi-major axis has not bee n mea - 
sured directly, but can be constrained. IWolff et all (|2002T) 
find eclipse durations 497.5 ± 6 s, which gives a eclipse half- 
angle 6 X = (6.543 ± 0.015)°. This constrains the inclination 
via Equation ((9j). Roche lobe overflow (Equation [7} implies 
a mean density for the donor star of 7.5 g/cm 3 , suggesting a 
main sequence or evolved companion. Maximum donor mass 
corresponds to a main sequence donor, with Md = O.42M0. 
Donor m ass can be lower if the companion is evolved, so 
following iHvnes et af] £2006) we take a minimum plausible 
companion mass of 0.07 Mq. For the NS we consider masses 
in the range 1.2 — 2AMq. This implies that a x sini lies in 
the range 0.11-0.84 lt-s. 



3.5 Kilohertz QPO sources 

These remain the most difficult sources, because of the un- 
certainty in the precise relationship between kHz QPO sep- 
aration (which varies) and spin frequency. Sources where a 
wide range of accretion rates have been sampled show varia- 
tion, those where only a few accretion rates have been sam- 
pled (including the pulsars that have kHz QPOs) do not. To 
gauge the uncertainty, consider the kHz QPO separations 
recorded for those sources where we have either a spin fre- 
quency or a burst oscillation frequency, illustrated in Figure 

El 

• Aql X-l: 550 Hz i ntermittent pulsar, twin kHz QP O 
separation 278 ± 18 Hz jBarret. Boutelier fc MUle3l2008h . 

• SAX J1808. 4-3658: 401 Hz pulsar. Twin kHz QPOs 
observe d once during the 2002 outburst. IWiinands et al.l 
ll2003al) report a separation of 195 ± 6 Hz, late r analysis by 
Ivan Straaten. van der Klis fc Wiinandsl (|2005l ) using a dif- 
ferent fitting technique gives a separation of 182 ± 8 Hz. 

• XTE J1807-294: 191 Hz pulsar. Twin kHz QPOs de- 
tected on several occasions during the 2003 outburst. The 
measured separations were consistent with the spin fre- 
quency, although the weighted averag e separation of 205 ± 6 
Hz exceeds the spin at the 2.3<r level (jLinares et alj|2005h . 

• 4U 1608-522: 619 Hz burst oscillations. Twin kHz 
QPOs observed in the 1996 and 1998 outbursts, separations 
in the range 225 - 325 Hz (|Mendez et al.lll998al fbTl. 

• SAX J1750. 8-2980: 601 Hz burst oscillations. Tenta- 
tive dete£t«3nj^wm_kHz QPOs with a separation of 317±9 
Hz (|Kaaret et alj|2002h . 

• 4U 1636-536: 582 Hz burst oscillations. Twin 
kHz QPOs detected on multiple occasi ons, with sep- 
aratio ns v arying from 240 — 325 Hz ( Wiinands et al.l 
1997al: 



Mendez. van der Klis fc van Paradiisl 



Jonker et alj 



2002 



Barret. Olive fc Millerll2005h 



di Salvo. Mendez fc van der Klia 



199S 
2003 



• KS 1731-260: 524 Hz burst oscillations. Twin kHz 
QPOs with separation 260. 3± 9. 6 Hz seen in one observation 
jWiinands fc van der Kli3ll99Tt ). 

• 4U 0614+09: 415 Hz burst oscillations. Twin kHz 
QPOs have been seen on several occasions (|Ford et al.l 
1 19971 ; Ivan Straaten et"al]|200Cl . l2002t ). Separations vary from 
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Figure 3. A comparison of twin kHz QPO separation and spin 
frequency or burst oscillation frequency for those sources that 
show both phenomena. The dotted line indicates equality of the 
two measures. For some objects kHz QPO separation is consistent 
with being constant: these are shown as single points with error 
bars. Note however that this may be due to poor sampling of 
source states. For five objects kHz QPO separation varies: these 
are shown as two points with a line indicating the range. 



238 ± 7 Hz up to 382 ± 7 Hz (perhaps even as low as 213 ± 9 
Hz although this figure is tentative). 

• 4U 1728-34: 363 Hz burst oscillations. Twin 
kHz QPOs detected on multiple occ asions, with sep- 
arations in the range 275-350 H z (IStrohmaver et al.l 



19961: iMendez fc van der Klisl Il999l : Idi Salvo et al.1 l200ll : 
Migliari et al.ll200St ). 

• 4U 1702-429: 330 Hz burst oscillations. Twin kHz 
QP Os detected once in 1997, with a separat ion of 333 ± 5 
Hz jMarkwardt. Strohmaver fc Swanklll999al ). 

• IGR J17191-2821: 294 Hz burst oscillations. Twi n 
kHz QPOs with separation 330 Hz (|Klein-Wolt et al.ll2007ch . 

• 4U 1916-053: 270 Hz burst oscillations. Twin kHz 
QPOs detected several times in 1996. Separation was con- 
sistent with being constant at 348 ± 12 Hz on four occasions; 
on the fifth it was 289 ± 5 Hz (|Boirin et al.ll2000l ). 

It is far from clear that there is a direct (or indeed any) 
relationship between kHz QPO separation and spin, partic - 
ular ly for the high frequency sources. See lYin et all (|2007m 
and lMendez fc Bellonil (|2007l ) for an extended discussion of 
this issue. 

We now summarize the properties for the kHz QPO 
sources. Note that we exclude from our data set the 
peculiar X-ray binary Cir X-l (1516-559), since it is 
not clear that our accretion torque m odel applies. Al- 
thou gh this source has twin kHz QPOs (|Boutloukos et all 
I2006T 1. it is thought to be a high mass X-ray bi- 
nary (jjonker. Nelemans fc Bassal 20071) with a highly 

>so|; lc 



eccen tric orbit (jMurdin et al.l 1 1980l ; lOosterbroek et all 
1 19951) where accretion disk formation is only sp oradic 



found that the long-term average flux is 1.1 x 10 ergs 
cm -2 s _1 . The best position for this source, from the opti- 
cal counterpart, is RA = 21 h 42 m 36'!9 1, Dec = +38°05 '27'. / 



( Johnston. Fender fc Wull 19991 ; I Johnston et al.ll200ll ). 



with an accuracy of_0"5 (B1950) (|Giacconi et al.l 1 19671 : 
iBradt fc McClintoclj[l983l). No more up to da te high preci- 
sion position is available. IGallowav et al] (|2008l ) use the PRE 
bursts to estimate a distance of 10 ± 2 kpc or 14 ± 3 kpc de- 
pending on composition. However optica l observations sug- 
gest a distance of onl y 7.2 ± 1.1 kpc (|Orosz fc Kuulkera 
1 19991 ). IWiinands et al.1 (|l998al ) detected twin kHz QPOs 
in RXTE data from a few hours on July 2 1997, when 
the source was on the horizontal branc h. Peak separa- 
tion was 346 ± 29 Hz. iKuznetsovl (|2002l ) re-analysed the 
same data and found a separation of 366 ±18 Hz. Un- 
successful searches for persistent pulsa tions were carried 
out using Ginga da t a from 1987-1989 (I Wood et all Il99ll : 
IVaughan et~aHll994! ). lKuulkers et all l| 19951 ) set upper lim- 
its on the presence of burst oscillations for frequenc ies be- 
low 25 6 Hz using EXOSAT data from 1983-1985, and lSmalei 
(| 19981 ) set upper limits on burst oscillations in the 200- 
600 Hz range for one burst observed with RXTE. The bi- 
nary orbit can be constrained by optical observations of 
the companion star V 1341 Cyg. The most recent spectro - 
scopic observations, bv lCasares. Charles fc Kuulkersl (|l998l ) 
give P orb = 9.8444 ± 0.0003 days, T 90 = 2449339.50 ± 0.03 
HJD, and e = 0.024 ± 0.015 (consistent at the 2a level with 
being zero, which we assume in our initial analysis). Error 
bars are lcr uncertainties. By measuring the projected veloc- 
ity of the secondary star and its rotational broadening, then 
assuming tidal locking and Roche lobe overflow, they infer 
D x sini = 29.9 ± 3.6 km/s. We use this value in our initial 
assessment of detectability. 



3.5.2 GX 340+0 (1642-455) 

This is a persistent Z source that has not shown X -ray 
bursts. Based on spectral modelling by I Ford et al.l (|2000l ) we 
set 25 ASM cts/s = 2.3 x 10~ 8 ergs cm" 2 s" 1 (bolometric) . 
The long-term average countrate is 30 cts/s. The best posi- 
tion, given by the radio counterpart, is RA = 16 h 45 m 44°60± 
0:02, Dec = -45°37 / 5 3"6± / .'3 (J2000) jPenninx et al.lll993l ). 
IChristian fc Swan 

3 (|1997| ) use Einstein observations to es- 
tablish an upper limit on the distance of 11 ± 3.3 kpJH. 
Twin kHz QPOs have been detected in RXTE data from 
199 7-8, with separations in th e range 275±24 Hz to 413±21 
Hz jjonker et al.lll998l . l2000r i. However, for the sample as a 
whole, the separation is formally consistent with being con- 
stant at 339±8 Hz. Unsuccessful searches for p ersistent pul- 
sations have been carried o ut using Ginga data (|Wood et al.l 
Il99ll : FVaughan et al.lll994l ). The orbital parameters are not 
known. 



3.5.3 4U 1735-44 

This is a bright persistent atoll source that has both 
bursts and superburs ts. Using the spectral modelling of 
IGallowav et~al] (|2008l ) we set 18 ASM cts/s = 9.1 x 10 -9 



3.5.1 Cyg X-2 (2142+380) 

This is a per s istent Z source that has X-ray bursts. 
IGallowav et al.l (|2008l ) carried out spectral modelling and 



7 iFord et al. I J_'t II M.| 1 quote a lower distance for this source, but 
the value given does not tally with that in the original reference 
that they cite. 
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ergs cm 2 s 1 (bolometric). The long-term average coun- 
trate is 14 cts/s. The best position, for the optical counter- 
part V926 Sco, is RA = 17 h 35 m 19!28 Dec = -44° 25' 20'.' 3, 
uncertainty < l" (|Bradt fc McCIinto"cklll983l ). No more up 
to date position is available. The distance inferred from the 
radi us expansion bursts ( assuming He limit) is 8.5 ± 1.3 
kpc (jGallowav et al.ll200sh . Twin kHz QPOs were seen m 
RXTE observations on May 30-31 1998. Peak separation 
varied from 296 ± 12 Hz up to 341 ± 7 Hz, inconsisten t 
at the 3.1a level with being constant (|Ford et alj Il99ct ) . 
iJongert fc van der Klisl { 1996) placed upper limits on the 
presence of burst oscillations (up to 256 Hz) using EXOSAT 
data: no figures have been reported for the RXTE burst sam- 
ple. The orbi tal parameters can b e constrained by optical 
observations. ICasares et al. I (|2006h report a recent spectro- 
scopic ephemeris, with P or b — 0.19383351 ± 3.2 x 10~ 7 days, 
in agreement with the most recent photometric ephemeris 
jAugusteiin et al1ll998l ). The time of inferior conjunction 
of the donor star T 90 = 24 5 281 3.495 ± 0.003 HJD. Further 
measurements and source modelling suggest that velocity 
semi-amplitude of the donor star is in the range 215-381 
km/s, with a mass ratio md/m x in the range 0.05-0.41. This 
would suggest u x sini < 156 km/s and we use this limit in 
our initial estimate of detectability. 



3.5.4 GX 5-1 (1758-250) 

This is a persistent non-bursting Z source, located in the 
highly absorb ed galactic bulge region. Based on spectral 
modelling by iFord et al.1 (|2000T l. we set 50 ASM cts/s = 
3.8 x 10~ 8 ergs cm~ 2 s _1 (bolometric). The long-term 
average ASM countrate is 70 cts/s. The best position, 
for the radio counterpart, is RA = 18 h 01 m 08!233, Dec 
= -25°04'42'.'044 (J2 000), positional uncertainty 0'.'040 
l|Berendsen et al.ll2000T ). Distance to this s ource is poorly 
constrained, although IChristian fc Swankl (|l997h give an 
upper limit of 9 ± 2.7 kpc0 Twin kHz QPOs have 
been reported in RXTE observations from 1996 - 2000 
(|Wiinands et al.1 Il998bl ; Ijonker et all I2002T ). Peak separa- 
tion is not constant, but varies from 232 ± 13 Hz up to 
344 ±12 Hz. Unsuccessful searches for persiste nt pulsations 
were carried out using Gin ga data from 1987 l|Wood et al.l 
Il99ll ; IVaughan et al.lll994T ). The orbital parameters are un- 
known. 



the globular cluster NGC 6624, and op tical observations im - 
ply a cluster distance of 7.6 ± 0.4 kpc jHeaslev et al.ll2000l ). 
The distance estimated from X-ra y bursts is ~ 6.6 kp c 
jVacca. Lewin fc van Paradiid Il986l ; iKuulkers et~ai1 120031 ). 
Twin kHz QPOs have been detected by RXTE on several 
occasions (ISmale Zhang fc Whitelll997t Izhang et al.lll998l ; 
lAltamirano et al" 20051 ). Measured separations vary in the 
range 220-350 Hz, although the uncertainties are such that 
the separation is co nsistent with being constant at « 275 
Hz. lDib et all l|2005h have placed extremely stringent limits 
on the presence of persistent pulsations up to 2000 Hz using 
archival RXTE dat a, for all a x s i n^ < 16.8 l t-ms (superced- 



ing ea rli er work by Wood et al.l (Il991 ) and IVaughan et al.l 
(| 19941 )). IJongert fc van der Klisl (|l996 ) placed upper limits 
on the presence of burst oscillations up to 256 Hz for EX- 
OSAT bursts. The most up to date X-ray ephemeris, us- 
ing all data from Ariel 5, Ginga and RXTE, gives P or b = 
685.0119 ± 1.02 x 10~ 4 s, P orb = (-7.54 ±3.21) x HT 13 s/s, 
and the time of maximum X-ray light is 2442803.63564 ± 
2.2 x 10" 4 HJD (|Chou fc Grindlay||200ll ). Noting that an 
ephe meris of similar accura cy can also be derived from UV 
data (|Anderson et al.lll997l ). we will assume that the refer- 
ence time could if required be related to a known phase in 
the orbit. There is also a super-orbital periodicity of ~ 176 
days, most likely due to perturbations of the orbital eccen- 
tricity (up to e = 0.004 by a third member of the system 
i|Chou fc Grindlay||200ll ; IZdziarski, Wen fc Gierhnskill2007l ) , 
but we neglect this in our initial assessment. The projected 
semi-major axis has not been measured, but can be con- 
strained. The assumption of Roche lobe overflow (Equa- 
tion [7} gives a mean donor density of 3300 g/cm 3 , im- 
plying a white dwarf companion. The donor must also be 
helium-rich to explain t he X-ray burst properties . Helium 
white dwarf models of iDelove fc Bildster] (|2003h suggest 
a donor mass in the range . 07 — 0.08Mp) (depending on 
temperature). [Anderson et al.l (|l997l ) inferred an inclination 
in the r ange 35 — 50° from obser v ations of the UV coun- 
terpart. iBallantvne fc Strohmaverl l|2004l ). analysing super- 
burst data, inferred a slightly lower value, so we will consider 
a minimum inclination of 30° . Then for neutron star masses 
in the range 1.2 — 2.4M© we p redict a x sin i in the range 7- 
20 lt- ms. Recent modelling by IZdziarski. Wen fc Gierliriskil 
(|2007l ) suggests that the super-orbital variability may even- 
tually pose even tighter constraints on the system. 



3.5.5 4U 1820-30 

This is a persistent bursting atoll source w ith regular dip- 
ping c ycles. Using the spectral modelling of lGallowav et all 
(|2008T ) we set 32 ASM cts/s = 2.1 x 10" 8 ergs cm" 2 
s" 1 (bolometric). The long-term average ASM countrate 
is 20.6 cts/s. The best position, from Hubble Space Tele- 
scope observations of the optical counterpart, is RA = 
18 h 23 m 40!4 53 ± 0^012, Dec = - 30°21'40'.'08 ± 0'.'15 (J2000, 
la errors), (|Sosin fc Kindll995lPI . The source is located in 



1989) 



2000) 



18 The distance of 6.4 ± 0.6 kpc de rived by iPenninx 
cited incorrectly as 9.2 ± 0.7 kpc by iFender fc Hendry 
iMigliari fc Fended ^0061) is highly model-dependent. 

19 Obtaining a more accurate radio position is com plicated by 
the presence of a nearby pulsar iMigliari et al.ll"2004l) . 



3.5.6 Sco X-l (1617-155) 

This the closest accreting neutron star in our study, and 
is a persisten t Z s ource. Using the spectral modelling 
of IFord et al l (|2000l ) we set 920 ASM cts/s = 4 x 10" 7 
ergs cm -2 s _1 (bolometric). The long-term average ASM 
countrate is 892 cts/s. The most accurate position for 
the source, given by VLBA measurements, is RA = 
16 h 1 9 m 55'!0850, Dec = -15°38'24'.'9, with an u ncertainty of 
0"5 (|Bradshaw. Fomalont fc Geldzahlerl Il999l ). The source 
has also has a measurab le proper motion, which we neglect 
in th is initial analysis (|Bradsha w. Fomalont fc Geldzahlerl 
Il999l). The distance, measured by paralla x, is 2.8 ± 0.3 kpc 
(|Bradshaw. Fomalont fc Geldzahler1ll999l ). Twin kHz QPOs 
are observed, with separ ati ons i n the range 240-310 Hz 
(Ivan der Klis et al.1 Il996l . Il997l ; iMendez fc van der Klisl 
|2000| ). There have been unsuccessful searches for 
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persistent pulsations up 



t o 25 6 
1986h_ 



Hz 
and 



using EXOSAT 
up to 512 Hz 



.,Middleditch fc Priedhorskv . 

using Ginga data |Wood et all Il99ll : iHertz et all 1 19921 : 

Photometric observations of the 
Sco imply an orbital pe- 



IVaughan et all ll994T ) 

optical counterpart V818 
riod of 68023.84 ± 0.08 
19751) 



although analysis 



s 

of 



[Gottlieb. Wright fc Lillerl 
R XTE ASM data by 
Vanderlinde. Levine fc Rappaportl (|2003h suggests that 
the true period cou ld in fact be sli g htly longer, at 



m 

68170 s. Analysis by ISteeghs" 

v x sin i = 40 ± 5 km/s. Assuming the orbital pe r iod o f 



Casares! (|2002h indicates 



iGottlieb. Wright fc Lillerl (|l975l ). ISteeghs fc CasareH (|2002 
derive an ephemeris with T 90 = 245 1 358.568 ± 0.003 HJE0- 



3.5.7 GX 17+2 (1813-140) 

This is a persistent Z source with X - ray b ursts. Using the 
spectral modelling of iGallowav et all <|2008l l we set 40 ASM 
cts/s = 1.6 x 10 -8 ergs cm -2 s _1 (bolometric). The long- 
term average ASM countrate is 45 cts/s. The most accurate 
position, from VLA observations of the radio counterpart, 
is RA = 18 h 16 m 1=389 ± 0'.'06, Dec = -14° 02' 10'/ 62 ± 0"0 4 
(J2000), lcr uncertainties in position (jPeutsch et al.|[l999h . 
This i s within the Chandra error circle for the X-ray po- 
sition (|Callanan et al.ll2002l ). Analysis of X-ray bursts sug- 
gests a distance of ei ther 9.8 ± 0.4 kpc or 12.8 ± 0.6 kpc 
(Galloway et alj|2008h . However, there are questions over 
how to correct for the super-Eddingt on persistent flux, an d 
the true distance could be lower (|Kuulkers et al. I I2002D . 
Twin kHz QPOs have b e en observed on mult iple occasions 
(|Wiinands et all Il997bl : iHoman et all I2002T ). Separation, 
which varies from 239 ± 17 Hz up to 308 ± 14 Hz, is not con- 
stant at the 97% confidence level. Upper limits on the pres- 
ence of persi s tent p ulsat ions in Ginga da t a wer e reported by 
IWood et all (|l99ll ) and IVaughan et all (|l994r ). Upper lim- 
its on the presence of burst oscilla tions in EXOSAT and 
RX TE data have been re ported by iKuulkers et al.l (119971 ) 
and lKuulkers et al.l (|2002l ). The orbital parameters are not 
known. 



3.5.8 XTE J2123-058 

This is a transient bursting atoll source in the Galactic Halo, 
with one recorded outbur s t in 19 98. Using the spectral mod- 
elling of IGallowav" et al.l (|2008l ) we set 6.4 ASM cts/s = 
2.1 x 10~ 9 ergs cm -2 s _1 (bolometric). The average ASM 
countrate over the 50 day outburst was 2.7 cts/s, yielding 
a long-term average countrate over RXTE's lifetime of 0.03 
cts/s. The most accurate position, as measured by Chan- 
dra, is RA = 2 1 h 23 m 14°54, Dec = - 5°47'53'.'2 (J20 0, un - 
certainty 0'.'6) jTomsick et all 120041 ). iTomsick et all (|200ll ) 
infer a distance of 8.5 ±2.5 kpc, consistent with the distance 
of 9.6 ± 1.3 kpc inferred bv lCasares et all (|2002l ). Twin kHz 



QPOs were detected during one observation, with separa- 
tions in the range 2 55 ± 1 4 to 276 ± 9 Hz dHoman et all 
Il999l ; ITomsick et a,l.|[l999h . ICasares et al.1 (|2002l ) use spec- 
troscopic and photometric me asurements to d erive P or b = 
21447.6 ± 0.2 s (la errors). ITomsick et"aH (|2002l ). how- 
ever, derive P or b = 21442.3 ± 1.8 (lcr errors). This discrep- 
ancy, and resu lt s from earlier photometric measurements by 
ITomsick et""aH (|l999l ); IZurita et ail (|2000l ). have yet to be 
resolved. The most recent time of minimum optical light 
T 90 = 2451779.652 ± 0.0 01 HJD. The projec ted velocity 
v^srai is also constrained. ICasares et all (|2002f ) attempt to 
measure this directly , and find i> x smi = 140 ± 27 km/s. 
ITomsick et"aH (|200ll. |2002| ) measure the projected orbital 
velocity and rotational velocity of the companion. Assum- 
ing a Roche lobe filling and tidally locked companion they 
infer the mass ratio, and hence a projected orbital velocity 
for the neutron star of « x sin i = llOtag km/s. lShahbaz et al.l 
(|2003l ) use th e projected or b ital v elocity of the companion 
measured by ICasares et al.l (|2002l ) and use more sophisti- 
cated models of the system to give the mass ratio. The re- 
sulting projected orbital velocity is v x smi — 103ly 6 km/s 
(90% confidence). In our i nitial analysis , howe ver, we use 
the direct measurement bv lCasares et all (|2002l ). 



3.5.9 GX 349+2 (1702-363, Sco X-2) 

This persistent Z sou rce does not show X-ray bur s ts. Us - 
ing the modelling of IZhang. Strohmaver fc Swankl (| 1998h , 
we set 43 ASM cts/s = 1.4 x 10" 8 ergs cm" 2 s" 1 (2- 
10 keV). There is no detailed spectral modelling avail- 
able for this source, so we adopt a bolometric correc- 
tion factor of 2. The long-term average ASM countrate 
is 50.2 cts/s. The best position, from VLA measurements 
of the radio counterpart, is RA = 17 h 02 m 22!93, Dec = - 
36°21 / 20 / /3 (B1950, accuracy C'5) l|Cooke fc Ponmanlll99ll ). 
Ilaria et al.1 (|2004l ) infer a distance of 3.6-4.4 kpc from 
BeppoSAX observations. Twin kHz QPOs were detected 
by RX TE in January 1998, with a separ ation of 266 ± 
13 Hz JZhang. Strohmaver fc Swankl Il998h . confirmed by 
lO'Neill et al.l(|2002l ). The source has however rarely been ob- 
served in the state where kHz QPOs are prevalent. A search 
for pers istent pulsations in G inga data from 1989 was unsuc- 
cessful (|Vaughan et al.|[l994l ). The binary period, measured 
using optical photomet ry and spectroscopy, has been the 
subject of some debate (ISouthwell. Casares fc Charleslll996l ; 
IWachter fc Margonlll996l ; lBarziv et a.l.lll997l ^) 7 but is now es- 
tablished as P or b = 22.5 ±0.1 hours (la error) (|Wachterl 
Il997l ). The other binary parameters have not been mea- 
sured, although we can constrain a x sini. The assumption 
of Roche lobe overflow (Equation [7} gives a mean donor 
density 0.2 g/cm 3 , which requires a donor that has evolved 
off the main sequence. In the absence of better constraints 
we will assume a mass ratio md/m x < 0.8. For the range of 
neutron star masses considered (1.2 — 2.4Mq) this implies 
a x sin i < 7 lt-s. 



20 lAbbott et al.l ll2007al) use a slightly larger uncertainty on Tqq 
in their searches for gravitati onal waves from Sco X-l, to account 
for the time elapsed since the lSteeghs fc Casares] j2002T) measure- 
ment. For consistency with the rest of our analysis we use the 
smaller uncertainty, assuming that the measurement could be re- 
done today to the same level of accuracy. 



4 SEARCHING WITH PARAMETER 
UNCERTAINTIES 

Let us summarize briefly the gravitational wave emission 
from a neutron star in a binary system. We assume that 
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the centre of mass of the binary is not accelerating in the 
solar system barycentre (SSB) frame. The timing model for 
the arrival times of the wave-fronts of the GW is taken to 
be the usual one jTavlor fc Weisberd Il989h . Let T be the 
arrival time of the wave at the SSB, r the proper time of 
emission in the rest frame of the neutron star, and t the 
time in the rest frame of the gravitational wave detector. 
The quantities T and r are related by 

T - T = t + A R + As + As (10) 

where A_r is the Roemer time delay accounting for the light 
travel time across the binary, Ae and As are respectively 
the orbital Einstein and Shapiro time delays in the binary, 
and To is a reference time. There are no additional timing 
delays due to dispersion. It turns out that, for our purposes, 
the Roemer delay is the most significant contribution. If r 
is the vector joining the centre of mass of the binary system 
with the neutron star, and n is the unit vector pointing from 
the SSB to the source, then 



There is then a similar relation between T and the arrival 
time t at the earth based detector, and we assume that this 
can be corrected for since the sky-position is known. 

In the models that we are considering, the intrinsic grav- 
itational wave frequency v depends on the spin frequency v 3 . 
The phase of the GW at the SSB is 

<t>{t) = $ + (12) 

where 

$(T) = 2-kv (t - - ' ^ \ . (13) 

Inclusion of frequency derivatives in this phase model is 
straightforward, and we do not write it down explicitly. Since 
the gravitational wave amplitudes are expected to be very 
weak and the output of the GW detectors dominated by 
noise, knowledge of the waveform, especially its phase, is 
crucial for detection. 

The phase Q(t) depends on the orbital parameters in- 
troduced in Section [3.21 P or h, Porb, Tasc, a x sini, and e. In 
addition, there are 2 parameters specifying the orientation 
of the orbital plane: the inclination angle i and the argument 
of periapsis uj. Of these 7 parameters, only 6 are required 
to define the phase mo del because of the projection a long 
the line of sight n; see iDhurandhar fc Vecchiol i|200lf ) for 
further details. Taking the spin frequency v s and its time 
derivative v s into account, we therefore have a total of 8 
parameters which determine the frequency evolution of the 
signal: (v 3 , i> a , a x sini, e, P or b, Porb, Tasc, w). 

This is clearly a very large parameter space, and a 
search over all these parameters using a sufficiently large 
data volume will be a big data analysis challenge. Let 
us therefore make some simplifying assumptions: u s — 0, 
Porb — 0, and e = 0; i.e. we assume a neutron star spin per- 
fectly balanced between accretion and gravitational radia- 
tion, and a circular orbit which does not decay appreciably 
over the course of the observation time. These assumptions 
may not hold for the sources and for the large observation 
times that we are considering, and an actual search might 
very well have to take some or all of these effects into ac- 
count. However, for assessing the detection prospects as we 



want to do here, this simplification is useful, since adding 
the extra parameters will further increase the number of 
templates. For some of the more promising sources at or 
near the detection threshold, these assumptions will need to 
be revisited in greater detail. Some of these extra parame- 
ters may need to included, and the resulting search might 
again become computationally difficult; this will be studied 
in greater detail in future work. 

In the case when the orbit is circular (e = 0), which 
we shall assume in the rest of this paper, the argument of 
periapsis and the initial orbital phase combine additively 
into a single parameter so that we are left with only 4 search 
parameters: A — (u s , a x sini, P or b, T asc ); we shall denote the 
components of A by A 8 with i = . . . 3. 

4.1 Template counting 

To determine the computational cost involved in searching 
the parameter space described above, we need to calculate 
the number of templates required. A calculation of the re- 
quired number of templates to search a portion of the pa- 
rameter space is based on demanding a certain maximum 
mismatch between the templates at neighbouring points in 
parameter space. This also guarantees that the true signal 
will not have more that the given mismatch to at least one 
of the search templates. The mismatch between waveform 
templates is measured simply as the fractional loss in the 
signal-to-noise ratio (SNR) when one waveform is filtered 
(folded) by the other. This fractional loss can be regarded as 
a distance measure between points in parameter space, and 
this leads naturally to the definition of a parameter-space 
metr i c gg I Sathvap rakash fc Dhurandharlll99ll . Il994l ; lOwenl 
Il996l : [Prrj3l2007bf l . Using the metric, we write the proper dis- 
tance squared (the "mismatch" ) between two infinitesimally 
separated parameter space points as 

m = ;/. '/A i7A (14) 

The size of the parameter space is then given by the volume 
measure determined by the metric in the usual way. Then, 
assuming that we cover this parameter space by a lattice 
grid, the number of templates is the total parameter space 
volume divided by the volume of each unit cell which makes 
up the lattice. 

The optimal choice of the lattice is determi ned by a so - 
lution to the so called sphere-covering problem (|PrixH2007r i. 
For our purposes, we shall use a simple cubic grid, and there 
are two reasons why it is acceptable to use this approxi- 
mation. First, the dimensionality of the reduced parameter 
space that we are looking at is low enough that the im- 
provement in the template pl acement ef ficiency is not more 
than a factor of about 2 or 3 <|Prixll2007l ). Furthermore, this 
improvement in the efficiency does not actually lead to a 
corresponding factor of 2-3 improvement in the sensitivity; 
the gain in the coherent integration time afforded by this im- 
provement is much smaller because the computational cost 
typically scales as a large power of the coherent integration 
time. Finally, the size of each unit cell is chosen based on 
the fractional loss in SNR, i.e. the mismatch m, that we are 
willing to tolerate; we shall use a reference value of m = 30% 
in this paper. 

The first detailed study of the parameter space met- 
ric for a neutron star in a binary orbit was carried out 
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in iDhurandhar fc Vecchid (|200l[ ) for a coherent matched 
filter search. The sea rch for gravit a tional radiation from 
Sco X-l reported in lAbbott et al.1 (|2007al ) was the first 
and so far, only application of this study. The aim of this 
section is mainly to collect some template counting equa- 
tions for later use. These equations can all be derived in 
a more or less straight f orwar d manner from the results of 
Dhurandhar fc Vecchio (12001 ) . The main difference is that 
Dhurandhar fc Vecchid 1 20011 ) use notation and variables 
targeted towards gravitational wave data analysis, while 
here we choose to use notation more familiar to an astron- 
omy/astrophysics audience. 

The first issue is the number of parameters which must 
be searched. Let us denote by AA 8 = A l nax — A^in the uncer- 
tainty in A 1 from astronomical observations; we assume the 
region to be rectangular. Since the proper length of the line 

in the A 1 direction is J^; max ^/giid\ % , a useful upper bound 
on this proper length is to use the maximum value of gu, i.e. 
(ffii ax ) 1//2 AA\ This proper distance can then be compared 
with our reference mismatch m, and we get a measure of the 
number of templates required in the X 1 direction: 



AA ! 



(15) 



If N X i < 1, it indicates the A 8 direction can be covered by 
just a single template, and the effective dimensionality of our 
parameter space is reduced by 1. The number of templates 
for the frequency are: 



N v = Av 



kT od 



(16) 



The uncertainty Av relates directly to the uncertainty in 
spin Av s . We shall ignore the correlations of v with the other 
parameters; this approximation will suffice for our purposes. 
For the other directions, the expressions for N X i have simple 
expressions in two regimes: T obs <C -Porb and T obs 3> P or b- 

In the limit of large observation times, Tobs Porb, we 
have: 



No. 



Nf 



n/2- 



AT a; 



AP G1 



im 



Dmin 
-'orb 



/o„2 f) max 1 
V07T a p i 



Tobs 



In the other limiting case T b s <fC Porb, we get 



N a 



Nj 



N P 



= Aa r 



= AT„. 



= AP 01 



2lV ^maxT^i 



yi75^(P- b n ) 4 
47r 3 a p ^ max To bs 



(17) 
(18) 
(19) 

(20) 
(21) 
(22) 



For a search over all three parameters at once, the total 
number of templates in (a p , P orb , Tasc) space is, for T obB 3> 

-Porb • 



A/ a 



L A[a p ]A[P^ 



A[T asc ] 



(23) 



while for T obs <C P orb it is: 



a P P orb T asc 



max ob | A[a P ]A[P o ^ ]A[T asc ] 
V702761 1500m 3 



(24) 

Note that iV ap iVp orb N Tasc # N ap p olb T asc • This happens be- 
cause the correlations between the different parameters (i.e. 
the off-diagonal terms in the metric) can be very important, 
especially for short observation times. 

Similarly, we shall require the equations in the case 
when one of the coordinates can be ignored and the search 
can be performed in a 2-dimensional subspace. The equa- 
tions for T ohs > Porb are: 



3 2 
^max 



A{a 2 p ]A[T as 

PorbTTl 

,~4/ max\2 2 rp 

k (a p ) f m „iobs 



A[P C 



orbJ A [Tasc] , 



'-A[al]A[p- D } 



(25) 
(26) 
(27) 



Finally, for T obs < P orb we get 



c\ 8 2 rp5 

— -A[a p ]A[T asc ] , 
45V35mP o b rb 



4tt v max i obs (a p J Arn-6 



135v 35m 

.2 rji6 



Na p P„ b 



8tt f max T obs rD- 6 i a r 2 

- A [-Porb] A K 



2835^5m 



(28) 

■Pasc] , (29) 

(30) 



We take the total number of templates to be A^ (if it 
exceeds unity) multiplied by the number of templates in 
(a p , Porb, Tasc) space. 

While these equations might not seem very illuminat- 
ing, two important features are worth remembering. First, 
and probably most importantly, the scaling of the number 
of templates with T obs is very different in the two regimes 
Pobs -C Porb and T obs 2> P or b- For example, in Equations 
([23} and ([24]), the scaling is C(T obs ) when T obs > P or b, 
while it is 0(To ba ) for small T obs . This will have important 
consequences for GW data analysis, as we shall see later. 
Second, computational cost issues become more important 
at higher frequencies and for tighter orbits, because the num- 
ber of templates typically increases when v m „ increases or 
Porb decreases. 

Finally, we say a few words about the positional ac- 
curacy required for the GW searches. We do not wish to 
consider searches over sky-position, and thus it is impor- 
tant to know the position sufficiently accurately before- 
hand. The periodic wave searches get their sky-position in- 
formation from the Doppler pattern of the frequency evo- 
lution. The sky-position accuracy A8 depends strongly on 
the coherent observation time T cob . For short observation 
times (< lyr) A6 increases as roughly 0(T o 2 bs ) or 0(T 3 ba ) 
(|Bradv et al.lll998l : |Prixll2007bl ). This increase in the sky- 
resolution eventually saturates when T co h becomes compa- 
rable to a year. The ultimate limit on A9 is the diffrac- 
tion limit, with Earth's orbit being the aperture size. Thus, 
the smallest error box for position will be A9 ~ A gw /1AU 
where A gw is the wavelength of the GW. This corresponds 
to about 8" at 50 Hz and it is inversely proportional to fre- 
quency; this requirement is easily met for all the sources 
we are considering. It might in fact be possible to use a sky- 
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position mismatch as a veto to rule out potential candidates. 
For a candidate with given values of the frequency and or- 
bital parameters, we could calculate the detection statistic 
at mismatched sky-positions and ve rify that the SNR does 
decrease as expected . The work of (Prix fc Itohll2005h (see 
also Appendix A of iKrishnan et al. ( 2004 )) which studies 
correlations in frequency and sky-position mismatch might 
be useful for this purpose. 



5 FUTURE DETECTABILITY 

5.1 Best case detectability for various emission 
models 

We will start by looking at the best case scenarios for de- 
tection, and ask what would happen if we knew all of the 
parameters to sufficient accuracy that we only had to search 
one template for each source, using the flux information col- 
lected in Section [3] and summarized in Table [1] In this case 
computational cost is not an issue, and we can integrate 
for long periods. We will assume that the system is in per- 
fect spin balance (so that we can neglect spin derivatives), 
and that gravitational wave torques are the only negative 
torques operating in the system. This means that we ne- 
glect any possible spin-down effects due to the interaction 
of the neutron star magnetosphere with the accretion disk, 
or any magnetic dipole spin-down. 

In Section [2] we gave an overview of the spin balance 
model and calculated the best case detectability assuming 
gravitational wave emission due to a 'mountain', balancing 
the long-term average flux. The results were shown in Figure 
[2] We also however need to consider whether the quadrupole 
Q required for spin balance is feasible. 



7.4 x 10 35 g cm 2 



R 3 \ 1/4 



Mi, 



rl F 1 / 2 
a kpc^ _g 



1 kHz 



5/2 



(31) 

where dk pc = d/(l kpc). Figure 2] shows the required 
quadrupole (scaled by 10 45 g cm 2 , the approximate mo- 
ment of inertia of a neutron star), assuming the long- 
term average flux. We also show the maximum sus- 
tainable quadrupole for an accre ted crust computed by 
lHaskell. Jones fc Anderssonl (|2006h (note that these authors 
compute Q22 = (15/87r) 1,/2 Q). The values required for 
the slower spinning bursters are at the upper boundary 
of what is thought to be feasible, and if the spin of the 
kHz QPO sources is slower than the measured separations, 
the required quadrupole may exceed this value. However, 
magnetic confinement could support larger q uadrupoles 
(|Melatos fc Pavnell2005l ; iPavne fc Melatos! 120061 ) . 

It is clear from Figure [2] that if the kHz QPO sources 
do have spins in the range inferred from the kHz QPO sep- 
aration (or higher) then Sco X-l is the only source that is 
in principle marginally detectable by Enhanced LIGO. Pre- 
dicted amplitude would of course rise if the spin frequency 
were substantially lower, and if this were the case both Sco 
X-l and GX 5-1 could be within detectable range for En- 
hanced LIGO (although the inferred quadrupole would be 
large, approaching the maximum thought possible). 

For Advanced LIGO, several of the kHz QPO sources 
are in principle detectable for a single template search in 
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Figure 4. The quadrupole Q required for spin balance in the 
mountain scenario given the long-term average flux. Q is scaled 
by 10 45 g cm 2 (the approximate moment of inertia of a neutron 
star). The frequency at which the kHz QPO symbols appear is the 
centre of the measured range of separations. The solid line illus- 
trates how the requisite Q would vary for Sco X-l if this is not the 
spin frequency. The dashed line show s the maximum feasible non- 
magne tic quadrupole calculated by lHaskell. Jones fc Anderssonl 
ll2006t> . The uncertainty in the inferred quadrupole is not shown, 
but can be substantial since it depends on the distance to the star 
(see Table [TJ. 



the broad band configuration. Two of the burst oscillation 
sources (XB 1254-690 and 4U 1728-34) are also marginally 
detectable within the given narrow band envelope. Several 
other burst oscillation sources might be detectable if it were 
possible to push the narrow band envelope further down 
towards the thermal noise floor. If the proposed Einstein 
Telescope reaches its design specification, several more of the 
burst and kHz QPO sources might be reachable, although 
the pulsars remain undetectable. 

One of the major uncertainties in our modelling is 
the response to variations in accretion rate. In order to 
gauge this effect we therefore consider an alternative model 
for the transients, one in which the gravitational wave 
torque balances the accretion torque in outburst. This as- 
sumes that the gravitational wave emission mechanism re- 
sponds very rapidly to the accretion. Although response 
timescales are not well studied, this scenario is not unrea- 
sonable - an accretion-induced mountain or unstable oscil- 
lation, for example, may well grow during outburst and de- 
cay during quiescence. Figure [S] compares the best case de- 
tectability for the transients if we balance the outburst flux 
rather than the long-term average flux (Table [1} ■ The in- 
tegration time in outburst is taken to be either the typ- 
ical outburst duration, or the maximum integration time, 
whichev er is larger. The required quadrup oles are all be- 
low the lHaskell. Jones fc Anderssonl (|20 06!l limit. Many of 
the sources show an improvement: the reduction in T b a 
is more than compensated for by the increase in flux. The 
burst oscillation sources X 1658-298 and KS 1731-260, for 
example, which are undetectable if we consider time-average 
flux, lie within this scenario on or just above the Ad- 
vanced LIGO Narrow Band envelope. The intermittent pul- 
sar HETE J1900. 1-2455 also becomes a more promising 
source for the Eintein Telescope. Further theoretical con- 
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Figure 5. Best case detcctability (single template) for the transients, in the mountain scenario. The left-hand panels show the detectability 
if gravitational wave emission balances the long-term average accretion rate, for T ^ a = 2 years. The right-hand panels show the 
detectability if we if we consider the quadrupolc necessary to balance the accretion torque during outburst. T v, a for each source is now 
either the outburst duration or the maximum integration time (2 years), whichever is longer. In order to compare sources with different 
integration times, we show log 10 of the ratio of the predicted to the detectable amplitude. The three rows are for different detector 
configurations: Top - Advanced LIGO broad band; Middle - Advanced LIGO Narrow Band Envelope; Bottom - Einstein Telescope. 



sideration should clearly be given to the issue of torque re- 
sponse timescales. 

We also consider the situation that would result if 
the accretion torque is balanced by spin-down due to 
an internal r-mode rather than a mountain. We first 
make the assumption that the star has a spin rate and 
temperature such that it can sustain an unstable r-mode 
(so that it lies in the 'r-m ode instability window', see 
lAndersson fc Kokkotasl l200ll ). This assumption may not 
be warranted, particularly for some of the more slowly 
rotating stars. We also assume that the r-mode amplitude 
is steady: r-mode unstable stars may well experienc e duty 
cycle s with short-lived periods of strong spin-down (iLevinl 
199S ; lAndersson et aL I l2000l ; lHevlll2002l ; iKinnev fc Mendelll 
20031 '). but there are some sc enarios in which r-mode emis- 
sion ma y be pe rsis tent lAndersson . Jones fc Kokkotasl 



20021: IWagonerl 120021: iNwyar fc Owenl l200rj : 



Bondarescu. Teukolskv fc Wasserman 



20071 ). The grav- 
itational wave torque associated with the dominant 
I — m = 2 r-mode (I and m being the standard angular 
quantum numbers) is 



-3jn a aMR 2 /T g 



(32) 



l|Andersson fc Kokkotasl l200li). The frequency of the emit- 
ted gravitational radiation is not a harmonic of the spin 
frequency: for the I — m — 2 r-mode it is 4/3 times the spin 
frequenc-J"! The quantity a is the mode amplitude. For an 
n = 1 polytrope model of the neutron star, J = 1.635 x 10~ 2 
and the radiation reaction timescale r g is given by 



(33) 



(jAndersson fc K okkotas 2001). The associated luminosity is 
i5gw = NgwUm/m, where uj m = u sv is the angular frequency 



21 We assume that this calculation of the relationship between 
mode frequency and spin is exact. In practice there is some the- 
oretical uncertainty, which would further increase the number of 
templates to be searched. 
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Figure 7. The r-mode amplitude a required for spin balance 
in the r-mode scenario, using the long-term average flux for 
all sources (persistent and transient). The dashed line shows 
the m aximum (saturation) amplitude calculated by lArras et al,l 
1120031 ). The uncertainty in the inferred amplitude is not shown, 
but can be substantial since it depends on the distance to the 
star. The solid line indicates the r-mode amplitude that would 
be required for Sco X-l if the spin frequency differs from that 
inferred from kHz QPO separation. 
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In Figure[6]we compare the predicted and detectable ampli- 
tudes for the r-mode scenario for T b s = 2 years. 

We now consider whether the required r-mode ampli- 
tude is feasible. Assuming spin balance, the mode amplitude 
is 



a — 6.9 x 10" 



(MhRlo) 



"kpc^_g 
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(35) 



Figure [7] shows the amplitude required for spin balance for 
each source for the long-term average flux. The required am- 
plitudes should be compar ed to the maximu m saturation 
amplitude a s computed by I Arras et al.l |2003f ): 



a s a8x 10" 



VI kHz/ 



5/2 



(36) 



The values required for the two slowest spinning 
bursters exceed the amplitude limit, as would the kHz QPO 
sources if the spin is much slower than the measured separa- 
tions. Howe ver for all other so urces the required amplitude 
is below the I Arras et all (120031 ) limit. 

The gravitational wave frequency is lower within the r- 
mode scenario than in the mountain scenario. Spin balance 
therefore requires a higher GW amplitude, making the r- 
mode emission model rather more optimistic. For Advanced 
LIGO many of the kHz QPO sources lie above the broad 
band noise curve. Two of the burst oscillation sources with 
reasonable values of a, 4U 1636-536 and 4U 1728-34, sit well 
above the Advanced LIGO narrow band envelope, while 4U 
1702-429 and KS 1731-260 lie on the curve (XB 1254-690, 



which would also be marginal for the narrow band config- 
uration, has an unfeasibly high a ). For the Einstein Tele- 
scope, however, these sources are well above the noise curve, 
and the intermittent pulsar Aql X-l is also marginally de- 
tectable. 

If we consider the transient scenario (balancing out- 
burst torque rather than time-average torque) the r-mode 
model also makes rather more optimistic predictions than 
the mountain model. Three burst oscillation sources that are 
marginally detectable at best when balancing long-term flux 
now come within range of Advanced LIGO's Narrow Band 
configuration: XB 1658-298, KS 1731-260 and 4U 1608-522. 
By the time we reach the sensitivity of the Einstein Tele- 
scope the intermittent pulsar HETE J1900. 1-2455 may also 
be detectable in outburst. 



5.2 Data analysis techniques 

We consider two kinds of searches over our param- 
eter space (either of these methods can also be 
used as pa rts of an optimize d mu lti-stage hierarchi- 
cal scheme (iBradv fc Creightonl bOQOl ; ICutler et all 120051 ; 
iMeinshausen, Bickel fc Ricell2007r )): 

(i) A full coherent matched-filter search using all avail- 
able data, possibly from multiple detectors. This involves 
a demodulation of the data for a given parameter space 

pint using the so called ^-statistic ; details can be found i n 
Jaranowski. Krolak fc Schutzll998l ; ICutler fc Schutzl2005l ). 
Such a coherent demodulation statistic can be augmented 
by the technique of com bining sidebands as described in 
iMessenger fc Woanl (|2007l ) (a dapted from similar te chniques 
used in radio pulsar searches, iRansom" et all (|2003l )). 

(ii) A power folding method where the available data is 
broken up into smaller segments (often called "stacks" in 
the GW literature); the duration of each stack is T co h and 
Attacks denotes the number of stacks. If there were no gaps 
in the data, we would have T b s = A r s t a cksT co h, but duty 
factors of ~ 80% are more likely. Each segment is coher- 
ently demodulated using (i) and excess power from each 
segment is combined without maintaining phase coherence. 
Examples of such power fo lding methods are: a) stack-slide 
(|Bradv fc Creightonl I2OO0I ) where the relevant statistic is 
simply a weighted sum of ^-statistic values (at the appro- 
priate frequency bins) 



"stacks 



(37) 



or b) the Hough transform method (|Krishnan et al.ll2004T ) 
where one adds weighted binary number counts 



(38) 



where ntk) is either or 1 depending on whether the T 
exceeds a certain threshold. In each case, the weights u>k are 
chosen to optimize the sensitivity. 

For both (i) and (ii), the computational cost is propor- 
tional to the total number of templates iVtemp and to the 
amount of data available. Thus, for (i), it is approximately 
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Figure 6. Best case detectability (single template search) using the long-term average flux, for the r-mode scenario. Symbols and lines 
as for Figure [2] 



ATobsiVcoh for some constant A, while for (ii) it is propor- 
tional to iVstacksT co h-/Vtemp- The precise value of the sensi- 
tivity of these searches clearly depends on the details of the 
analysis method and software, and the quality of data, such 
as the duty cycle of the detector (which might reduce the 
amount of data actually available), the presence of noise 
artifacts such as spectral disturbances, uncertainties in the 
calibration of the detector, and so on. It is thus not possible 
to estimate the sensitivity to better than, say, ~ 5 — 10% 
without actually carrying out the search, and it is in fact 
even pointless to try and do so in this paper. It is of course 
possible to get semi-realistic estimates and this is what we 
shall do, but these uncertainties should always be kept in 
mind. 

Starting with the statistical factor, note that if the 
threshold corresponding to a single trial false alarm rate is a, 
the probability that the threshold is crossed at least once in 
■^trials independent trials is FA = 1 — (1 — q)^"^ 8 sa aN trMs 
when adVtriais <C 1 (|jaranowski. Krolak fc Schutd[l99sh . So 
we choose a — FA/NtriaXs to ensure that the total false 
alarm probability is FA. We equate the number of trials 
with the total number of templates that must be searched. 
Strictly speaking, this is not true because the different tem- 
plates are not completely independent. This is usually not 
a significant effect as long as the mismatch m used to con- 
struct the template bank is not too small; we shall use a 
reference value of in — 0.3. To make matters more compli- 
cated, while there are reliable calculations for the template 
counting for (i) as presented in Sec. 14.11 there are as yet no 



reliable estimates for (ii). We shall therefore consider only 
the number of templates in each coherent segment. Fortu- 
nately, the dependence of the statistical factor is not very 
steep (in fact slower than logarithmic) with the number of 
trials, so this does not make a significant difference to our 
results. 

Some details of the statistical calculation are in order. 
The choice of a determines the thresholds for both the co- 
herent jF-statistic search and the semi-coherent search. The 
exact relation is however different. The ^-statistic follows a 
^-distribution with 4 degrees of freedom (it is actually IT 
which is \ 2 distributed). For the semi-coherent statistic p is, 
to a reasonable approximation, simply Gaussian (assuming 
Attacks to be sufficiently large); it is actually \ 2 with 4A r sta cks 
degrees of freedom. The assumption of gaussianity may be 
questionable here because we are after all dealing with the 
tail of the distribution, and the central limit theorem may 
not be reliable especially when Attacks is not particularly 
large. A further complication arises when p is not simply a 
sum of the ^-statistic values but is perhaps a Hough statistic 
(|Krishnan et alj|2004h when the distribution of p is closer to 
a binomial (which can also be approximated by a Gaussian) . 
We shall nevertheless assume gaussianity for our purposes 
and leave a more detailed study for future work. The rela- 
tions between a and the corresponding threshold jT th and 
Pth are then different in the two cases: 

a = (1 + Fth)e-^ b , a = ierfc ( -^-) (39) 
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Figure 8. Dependence on the number of independent trials (tem- 
plates) of the statistical factor F s tat that appears in the de- 
tectability equation (1406 . Note that -Fstat = 11.4 when W tr ; a i s = 1 
which reproduces equation J6j. 



where erfc is the complementary error function, and a is 
the standard deviation of p. In the presence of a signal with 
amplitude ho, the mean values of T and p are increased by 
an amount proportional to the SNR 2 , and thus to h%. Tak- 
ing a fixed value of the false dismissal rate j3 determines the 
minimum value of ho needed to exceed the thresholds pth 
and jT th ; we shall always take (3 = 0.1 as the reference value. 
Folding in a uniform averaging over all possible pulsar orien- 
tations along with the statistical factor in a single parameter 
-Pltatj we get the following expression for the sensitivity 



/ill 



N. 



1/4 



S h 
DT co 



(40) 



where F s tat varies with the number of trials required, 
as shown in Figure [8] Th is factor was also d iscussed in 
(| Abbott et al.ll2007al 1 and (| Abbott et alj l2005al ) which are 
examples of (i) and (ii) respectively. In each case, Fstat in- 
creases slower than logarithmically with Aerials, and the sta- 
tistical factor is worse for the x 2 distribution as expected be- 
cause it is not as sharply peaked as the Gaussian. Finally, it 
is perhaps worth mentioning that averaging over the pulsar 
orientation cos t may not always be appropriate, especially 
if there happen to be independent estimates of the neutron 
star orientation. If the value of cos l happens to be anywhere 
in the top 1 — p percentile of a uniform distribution (i.e. if 
we average over all | cos t| ^ p), then 



F 3 t 



-^St a 



v/l + ll(p + p 2 )/16 + (p 3 +p 4 )/16 



(41) 



So for example, if p = 0.9, we get F sta t — > Fkat/1-5. This 
could be important for sources which happen to be near the 
detection threshold. 

So far we have incorporated the statistical and geo- 
metrical factors in Equation (|40|l . For template bank based 
searches, the sensitivity is further degraded because of the 
discreteness of the template grid. For a template bank cre- 
ated with a maximum mismatch of m, assuming a cubic 
grid, the average degradation in the SNR is (1 — m/3). This 



corresponds to a degradation of \/T 
0.95 for m = 0.3. 



m/3 in ho which is 



The next question we need to address is the computa- 
tional cost and whether it is necessary or worthwhile to do 
the semi-coherent search (ii). Clearly, for a given amount of 
data while we ideally want to make the coherent integration 
time T co h as large as possible, and the number of segments 
-^stacks correspondingly small, the choice of these parameters 
is dictated by the computational cost. This computational 
burden is in fact, by far, what limits the search sensitivity 
(As we just saw, Fstat has a very weak dependence on the 
number of templates; we need to change the number of tem- 
plates by orders of magnitude before it has an appreciable 
effect). The bigger effect of reducing the number of tem- 
plates by better astrophysical modelling or observations is 
that it allows us to have a larger coherent integration time 

T co h • 

For a given source, we would first like to know if a 
semi-coherent search would be useful. The obvious answer 
is: whenever computational cost is an issue, but a more 
quantitative answer is also easy to get. Let us compare the 
sensitivities of (i) and (ii) for a given computational cost. 
Let us assume that the cost for the coherent analysis scales 
as a power of T b s : C co h = AT„ bB for some constants A 
and k. So if we assume a fixed value of the computational 
cost Co, then T b s = (Co/A) 1 ^. Thus, the sensitivity is 
ho° h oc T^J 2 = (C /A)~ 1/2k . For a semi-coherent search 
on the other hand, we have Co ~ A Attacks T C oh (this is 
true if the cost of combining the different stacks is negli- 
gible). In this case we get T coh = {Co/AN st ^cks) lllk ■ Using 



I,- 

n 



oc AT 



i semi 



Lcoh 
n 



(2-fc)/4fc 



(42) 



This tells us that a semi-coherent search is not effective (i.e. 

^semi-coh > ^coh) for fc ^ ^ ^ j t getg mQre &nd mQre 

effective for larger k (this conclusion is robust: it is not af- 
fected by the approximation of neglecting the cost of the 
semi-coherent step). 

In our present case, recall from equations (|17p - (|28|l that 
we have very different scalings in the regimes T b s 3> -Porb 
and Tobs •C F or b. So, for binary systems for which the pa- 
rameters have been sufficiently constrained astronomically 
and we can afford large integration times, we shall assume 
that we are better off doing a simple coherent search. For 
others, a semi-coherent strategy has a much larger impact I 22 ! 
For each potential source, our strategy is to first estimate the 
number of templates that are required for a coherent search. 
With reasonable estimates of available computational re- 
sources, and assuming T b s = 2 years, this determines the 
maximum coherent integration time T co h that can be anal- 
ysed. If T co h is not much smaller than P rb, then by the 
previous argument we assume that the gain in using a semi- 
coherent search is not significant, and we restrict ourselves 
to a pure coherent search. On the other hand, if this T co h 
does turn out to be much less than P or b, then we consider a 
semi-coherent search and estimate T co h and Attacks (assum- 
ing here the cost of the semi-coherent step to be comparable 



22 Semi-coherent searches might start becoming more important 
as soon as we need to start including other parameters, especially 



frequency derivatives and P c 



orb ■ 
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Figure 9. The factor by which F s t a t increases over the single 
template value (-Fstat = 11.4) plotted against the predicted am- 
plitude for the mountain scenario (Equation[5]l. The brightest and 
most promising sources are the most poorly constrained, and this 
is reflected in the effect on -F s t a t • The difference in -F s tat between 
the two strongest sources, Sco X-l and GX 5-1, is primarily due 
to the fact that Sco X-l has much better orbital constraints. 

to the coherent step). Equation f (|40[) ) then yields an esti- 
mate of the search sensitivity hg cnB for each potential source. 
We take the statistical factor F a tat in Equation f (|40)l ) using 
just the number of templates in the coherent step. There 
is as yet no systematic study of the semi-coherent metric 
for binary systems, and this estimate should be updated as 
soon as these calculations are available. We expect that this 
approximation will not make a qualitative difference to our 
final results because of the weak dependence of F s t a t on the 
number of templates, but this needs to be verified. In the 
previous section, we estimated the estimated GW amplitude 
/iq XP from each source, and the ratio /iq XP /ho cns can be a use- 
ful "detectability" ranking for each source. 

In the cases when P or b is essentially unknown, the 
ranges of r asc and a x sin i that we have chosen depend on 
Porb ■ This means that the parameter space region is not rect- 
angular and thus the template counting equations of Section 
14. II will over-estimate the number of templates significantly, 
though not by orders of magnitude. 

5.3 Numbers of templates 

We can now compute the effect on our searches of param- 
eter uncertainty. We will start by considering the effect on 
the statistical factor F s tat , neglecting computational cost is- 
sues. In other words, we assume that it is feasible to do a 
coherent fold. In Tables 0E] and 3] we summarize the param- 
eter ranges assumed for the pulsars, bursters and kHz QPO 
sources respectively. The Tables show the associated N X i 
for each parameter (for T b s = 2 years), the resulting total 
number of templates that must be searched iVtemp, and the 
effect on -Fstat ■ Figure [9] shows the change in the statistical 
factor for each source. 

In Figure [2] were we showed results for a single template 
search, each source had F s tat = 11.4. This meant that we 
could plot a single detectability threshold curve for each de- 
tector. When each source has a different Fstat, however, the 
detectability threshold curves differ for each source. One way 



of comparing detectability for sources with different noise 
threshold curves is to plot the ratio of emitted to detectable 
amplitude: this is the approach that we took for the tran- 
sients (Figure [5}, and that we will adopt when we come to 
consider computational cost. What we can also do, however, 
to see the effect, is to scale the predicted amplitude ho by 
11.4/Fstat. Predicted amplitude does not really change, of 
course, but it is a useful way of visualising the impact of the 
statistical factor. 

Figure lTUI shows the impact on the mountain scenario: it 
should be compared to Figure [2] Only for the most tightly 
constrained pulsars is detectability largely unchanged: for 
the majority of sources the detectable amplitude falls by a 
factor of 2-3 compared to a single template search. Although 
this does not sound like a great deal, it is sufficient to push 
all sources except Sco X-l below the detection threshold for 
Advanced LIGO broad baneful. The effect on the r-mode sce- 
nario, summarized in Tables I2T41 is slightly less severe, leav- 
ing Sco X-l, GX 5-1, GX 349+2 and 4U 1820-30 above the 
detection threshold for Advanced LIGO broad band. The sit- 
uation is better for the narrow band configuration: although 
none of the burst oscillation sources remain within range (for 
any of the emission scenarios considered), several of the kHz 
QPO sources are still viable. The spins for these sources are 
poorly constrained. However, the anticipated narrow band 
configurations (see top panel of Figure [1]) have a reasonably 
broad bandwidth, leaving ample scope for searches. 

We can now take the final step and look at the impact 
of computational constraints on searches involving multi- 
ple templates. The available computational power sets the 
length of data T b s that can feasibly be analysed within a 
given amount of time. In Table[5]we summarize the impact in 
the situation where we assume a maximum analysis time of 
2 years, assuming that for Advanced LIGO we have compu- 
tational power 50 times greater than that currently available 
within the LSC, while for the Einstein Telescope we assume 
100 times more computing than at present. These assump- 
tions are, of course, arbitrary, since the computing power 
that can be applied to future searches depends not just on 
technology (e.g. Moore's Law) but also on improvements in 
data analysis techniques, and of course also on how much 
money is spent on computing by future projects. We have 
taken numbers we feel are defensible, but they may turn out 
to be conservative, especially for the Einstein Telescope. 

For the pulsars, computing constraints have little effect 
on the detectability of the sources; however, these are not 
likely to be detectable in any case. But for the bursters and 
the kHz QPO sources, computational constraints lead to a 
major reduction in the T bs that it is feasible to analyse. For 
many sources we have to resort to a semi-coherent search, 
as a coherent fold is no longer possible. Figures QT] and [12] 
show the effects on detectability for the mountain and r- 
mode scenarios respectively, taking into account both the 
change in F s t a t and the reduced observation times for our 
assumed computing power. 

The additional impact of computational limitation is 
substantial. For Advanced LIGO's broad band configura- 
tion, only Sco X-l and 4U 1820-30 (the latter in the r-mode 



Sco5 X-l would also be marginal for Enhanced LIGO, over a 
restricted frequency range. 
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Figure 10. Effect on detectability for the mountain scenario (assuming long-term average flux and a coherent fold with T t, s = 2 years), 
taking into account the effect on F s tat associated with the fact that TVtemp > 1. Compare to the best case detectability shown in Figure 
[2] As in Figure [2] the noise curves are computed assuming F s t a t = 11.4 (the single template value), but we have scaled the predicted 
amplitudes to reflect the fact that F s tat is larger. Although this is not strictly "correct" (it is the thresholds that should move, not the 
predicted amplitudes) this is a useful way to visualise the impact. See the text for more details. 



scenario) remain above the detectability threshold. These 
two sources are also the only two left above the narrow band 
envelope threshold, although there are other sources just be- 
low the envelope that might be detectable if the narrow band 
configuration were able to push closer to the thermal noise 
floor. 

Parameter uncertainty clearly poses a major problem, 
even for stars where we have some indication of the spin 
rat J^l. If the parameter space volume can be reduced, how- 
ever, then the statistical and computational restrictions will 
have less impact. In the following section we will examine 
this in more detail and consider how best to solve it. 



but better constrained burst oscillation sources. The best 
constrained sources, the accreting millisecond pulsars, are 
expected to emit at such a low level that they are unlikely 
to be detectable by second generation detectors. This is par- 
ticularly depressing since our calculations have been carried 
out within the context of a fairly optimistic modelling and 
analysis scenario. We have not as yet considered any addi- 
tional sources of spin-down, and have also neglected physi- 
cally reasonable complications such as spin variability which 
would increase the number of templates still further. Spin 
variability would of course also lead to decoherence of the 
signal for sources where we cannot track the spin, reducing 
integration times to the decoherence time Td CC oh, defined by 



6 DISCUSSION 

6.1 Current and future prospects 

It is clear from Section [5] that the detection of gravitational 
waves from accreting neutron stars is a difficult task. The 
X-ray bright kHz QPO sources suffer from parameter un- 
certainty, forcing us to look in addition at the much weaker 



This problem is of course even more pronounced for the ac- 
creting neutron stars where we have no indication of spin rate. 



Tdccoh^s = 1. (43) 

One can get some idea of the worst case scenario by calculat- 
ing the u 3 that would result if the source were spinning up at 
the maximal rate due to the accretion torque (Equation 1). 
The results are given in Table [5] for the bright kHz sources 
this worst case decoherence time could be ~ 1 week, and 
this is something that will need to be considered in future 
studies of detectability. 

In order to improve prospects for those sources that 
are in principle detectable, we clearly need to find ways of 
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Figure 11. Detectability for the mountain scenario (long-term average flux, maximum T t, s = 2 years), taking into account computational 
limitations as well as the effect on -F s tat ■ The computational cost of searching multiple templates affects the feasible integration time and 
the choice of coherent / semi-coherent search technique (Table O . In order to compare detectability for sources with different integration 
times we plot log 10 of the ratio of predicted to detectable amplitude. The three rows are for different detector configurations: Top - 
Advanced LIGO broad band; Middle - Advanced LIGO narrow band envelope; Bottom - Einstein Telescope. The left hand panels show 
detectability for the best case (single template) search, as in Figure [2] the right hand panels show the situation after including the 
statistical and computational limitation. For Advanced LIGO we assume a 50-fold increase in computational efficiency as compared to 
the current situation; for the Einstein Telescope we assume a 100-fold increase. In both cases we assume a maximum analysis time of 2 
years. 



improving source constraints and removing computational 
limitations. In the sections that follow we detail the actions 
that would lead to the biggest improvements. 

6.2 Astronomical observations 

The major obstacle to detection by the current and immi- 
nent generation of detectors is clearly parameter uncertainty, 
with spin uncertainty having the largest impact (see Ta- 
bles [2][4j) . Astronomical observations might help to constrain 
source properties. 

The single most valuable thing that could be done to 
improve the current situation is to determine more precisely 
the NS spin in the bright kHz QPO and burst sources. A 
substantial archive of high time resolution X-ray data exists 
for many of these sources from RXTE and its predecessors. 
However the only dee p search for pulsat ions in the literature 
is that carried out bv lDib et al] (|2005l ) for 4U 1820-30. Se- 
rious consideration should be given to carrying out similar 



analysis for all of the kHz QPO and burst sources, most par- 
ticularly Sco X-l. Thorough searc hes for intermittent p ulsa- 
tions, such as tho s e carr ied out bv lCasella et all (|2008l ) and 
lAltamirano et aL I (|2008f ). would also be worthwhile. 

It may also be possible to find burst oscillations in 
the kHz QPO sources: five of the sources that we have 
analysed are known to burst. If a sufficient archive of 
bursts can be built up from these sources burst oscilla- 
tions may well be detected either in individual bursts or 
by stacking power spectra from multiple bursts (the tech- 
niqu e used to find weak burst oscil lations from EXO 0748- 
676 (|Villarreal fc Strohmavei1l2004l ). We also need to verify 
the burst oscillation frequencies for the seven sources with 
tentative or single burst detections, in particular the four 
that are above the detection threshold for some of the sce- 
narios that we have examined: XB 1254-690, 4U 1916-053, 
MXB 1730-335 and 4U 0614+09. 

Identifying the orbital period can also make a substan- 
tial difference to the number of templates searched. Con- 
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Figure 12. As Figure ITU but for the r-mode scenario. 



sider for example the difference between the two kHz QPO 
sources Sco X-l and GX 5-1. Both are expected to be strong 
emitters, but Sco X-l suffers much less from parameter un- 
certainty than GX 5-1 because the former has a well con- 
strained orbital period. Three of the most promising kHz 
QPO sources have no orbital constraints (GX 340+0, GX 
5-1 and GX 17+2), while a fourth (GX 349+2) has only 
a relatively weak constraint. Several of the most promising 
burst oscillation sources also fall into this category: atten- 
tion should focus on 4U 1608-522, KS 1731-260, 4U 1728- 
34, 4U 1702-429, and MXB 1730-335, sources that could be 
detectable in some scenarios. Immense progress has been 
made in recent years in constraining orbital parameters for 
LMXBs using optical/IR observations, particularly with the 
Bowen technique of spectroscopy. Identification of counter- 
parts is often difficult, but there would be a substantial pay- 
off in terms of gravitational wave detection prospects. 

6.3 Astrophysical modelling issues 

Both burst oscillation frequency and kHz QPO frequency 
are proxy measures for the stellar spin rate. The precise link 
to spin for each measure is not clear because in neither case 
do we understand the mechanism. Astrophysical modelling 
to pin down the cause of the two phenomena could therefore 
shed light on the reliability and accuracy of the spin proxy. 



Several models have been suggested for burst oscilla- 
tions. These include global oscillations of the surface lay- 
ers of the neutron star, large-scale drifting vortices excited 
by the passage of the flame fronts, or hotspots linked in 
some way to persistent pulsations. All of the models have 
shortcomings, and progress towards understanding this phe- 
nomenon has to some degree stalled. The situation for the 
kHz QPOs is however even worse. Ther e are several mod- 
els, summarized in detail in Ivan der Klisl (|2006T ). All involve 
either orbital motions of material within the disk or disk os- 
cillations, and many also require some mechanism to select 
preferred radii within the disk. Developing testable predic- 
tions that would distinguish between the different models, 
however, is very difficult. A substantial body of data exists 
within the RXTE archive to test models if such predictions 
can be made. The identification of a robust link between kHz 
QPO frequency and stellar spin would have a major impact 
in reducing parameter uncertainty for the most promising 
sources for Advanced LIGO. 



6.4 Data analysis and detector issues 

In the absence of computational limitations, even taking 
into account the effect on statistics, there would be up to 
8 sources potentially detectable by Advanced LIGO (Figure 
UOp . We therefore need to improve the data analysis tools 
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used for the searches. We need significant improvements in 
all relevant aspects of the data analysis pipelines: acceler- 
ating existing software through improvements to both soft- 
ware algorithms and computer hardware, and also develop- 
ing other data analysis techniques. The software algorithms 
being used in the LIGO data analysis software have im- 
proved continually over the years since the first LSC publi- 
cation on periodic waves in 2003; we expect these improve- 
ments to continue. It is possible that we may get an improve- 
ment of more than a order of magnitude in existing codes 
over the next several years. The improvements in computer 
hardware, even just following Moore's law, will yield an ad- 
ditional factor of about 16 in the next 8 years, i.e. by the 
time we expect Advanced LIGO to be operating. Further- 
more, computing platforms like Einstein@Home allow us to 
increase the total number of computers available to do the 
analysis. On top of this, it might also be possible to use spe- 
cial purpose hardware for the periodic wave searches. Most 
of the analysis involves a large number of relatively simple 
operations, and it might be quite feasible to design chips 
which are efficient for these particular calculations. An ex- 
ample of these are the GRAPE special-purpose computers 
which have proven t o be extremely useful in astrophysical 
N-body simulations (|Makino fc Taiiilll995T ) and in molecu- 
lar dynamics calculations in condensed matter physics. It is 
hard to anticipate these developments, but an improvement 
of 4 or 5 orders of magnitude in effective computing power 
might be feasible. 

Apart from these technical improvements, it is also pos- 
sible to develop new analysis methods. A good example of 
this is the cross-correlation method used in the Sco X-l 
search (| Abbott et al.ll2007d ). This was neither a matched 
filter nor a standard time-frequency semi-coherent search. It 
was instead based on aperture synthesis, i.e. the fact that we 
have multiple detectors in operation, and that they all see es- 
sentially the same GW signal at any given time. This method 
has so far mostly been used in the stochastic background 
searches, and is now being adap ted to the periodic wave 
searches (|Dhurandhar et EiLl 12008). Apart from the compu- 
tational efficiency, this method is also relatively insensitive 
to the uncertainty in the signal model caused, e.g. by pul- 
sar glitches (though this leads to a correspondingly greater 
computational cost in any follow-ups that must be done for 
verifying a detection and for parameter estimation) . Further 
searches using this method, in combination with the other 
methods discussed in this paper seem to be quite promising, 
especially given that by the time of Advanced LIGO there 
will be a third interferometer of comparable sensitivity in 
Virgo. 

In the near term while these improvements are in 
progress, it seems reasonable to search for some of the 
brighter sources such as Sco X-l not at the frequency implied 
by the kHz QPOs, but rather at lower frequencies corre- 
sponding to the sweet spot of the detectors. Not only are the 
instruments more sensitive at these frequencies, but this is 
computationally easier and is astrophysically well motivated 
if the link between the kHz QPO and the spin frequency is 
not validated. 

We conclude this section with a brief discussion of pos- 
sible Advanced LIGO configurations and some of its impli- 
cations for our purposes. Figure [13] shows noise curves for 
three cases: an example broad band configuration with some 
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Figure 13. Some possible Advanced LIGO configurations. See 
text for discussion. 



tuning of the signal recycling cavity optimized for binary 
neutron star (BNS) inspirals, the envelope of possible nar- 
row band configurations, and the zero-detuned broad band 
configuration. The zero-detuned curve and the narrow band 
envelope are the same as in Fig.[T] The BNS curve is relevant 
because binary neutron star inspirals will be among the key 
targets for Advanced LIGO and the detector might possibly 
be operated in this configuration for significant durations. 
Thus, while BNS signals are not our concern in this paper, 
it is important to consider the impact of this configuration 
for our purposes. It is obvious that the BNS curve affects us 
adversely above ~ 500 Hz which, as we see from Figure [2] is 
precisely where the kHz QPO sources are expected to lie. 

These curves highlight the importance of spin measure- 
ments. For all potential sources with frequencies between, 
say, 150 and 2000 Hz, the narrow banding gains us no more 
than a factor of 2. The best case for narrow banding is then 
a compelling source whose frequency is very well known and 
is within this range and is not detectable by the broad-band 
curve. The other possibility is to employ the zero-detuned 
noise curve for detection and then use the narrow band curve 
as a follow up to confirm detection or to increase SNR for ex- 
tracting better astrophysical information. Given the present 
uncertainties in the spin frequency, this latter possibility 
seems to be the best option for detecting periodic waves 
from LMXBs. This conclusion could change if there turn 
out to be significant improvements in the narrow band noise 
curves. 

6.5 Conclusions 

We have shown that gravitational wave observations of ac- 
creting neutron stars will be challenging. It is therefore 
worthwhile considering what is to be gained from such an ef- 
fort: what new information gravitational wave observations 
will bring to the study of these sources. 

The input from electromagnetic observations that is re- 
quired for these searches is primarily on kinematical param- 
eters like spin rates and orbital ephemerides. The gravita- 
tional waves that are generated carry information on the 
orientation and dynamics of the neutron star and its inter- 
action with the accretion disk. The elliptical polarization 
of the gravitational waves, which will emerge from anal- 
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ysis of the detected signal, determines the orientation of 
the neutron-star spin axis. The degree of alignment of this 
with the inferred disk orientation, especially for the slower- 
spinning (presumably younger) systems, will constrain mod- 
els of binary evolution and the formation of the NS. 

In systems where observations and modelling constrain 
the NS spin to a relatively narrow range, the GW frequency 
will determine whether the r-mode or mountain scenario, or 
indeed some other scenario, is the appropriate one. The GW 
amplitude will then determine the degree of mass asymme- 
try (mountains) or the size of the velocity field (r- modes). 
In both cases this will open for the first time a wealth of 
opportunities for studying the physics of NS interiors. 

For all detected systems, the measured GW amplitude 
and frequency will tell us how much angular momentum is 
being carried away in GWs; by measuring or limiting the 
rate of change of the GW frequency we can then infer the 
rate at which angular momentum is being accreted from the 
disk. This will be an important constraint on models of the 
disk and the NS magnetic field. 

Clearly, if GW searches turn up signals in unexpected 
places in parameter space, this will challenge prevailing mod- 
els for these systems. For example, if the GW frequency 
equals the X-ray pulsation frequency instead of twice or four 
thirds its value, then this might imply that the pulsation fre- 
quency is actually twice the true spin frequency. 

Even in the worst case, where a sufficiently sensitive 
search fails to detect GWs at the expected amplitude, the 
negative result could eliminate the GW option for spin 
balance and demonstrate that this somehow resulted from 
a propeller-type torque from the NS's magnetic field. Of 
course, one would have to have confidence that the GWs 
were not coming out at a different frequency, so the search 
would have to include that as a parameter. 

The additional astronomical input that could make the 
difference between detecting or not detecting GWs could 
come from a variety of observations. Clearly, long-term X- 
ray timing is highly desirable, and here it is disappointing 
that the future of RXTE is limited and no X-ray mission 
with comparable capabilities seems to be planned for the 
period before Advanced LIGO comes online. Thorough ex- 
ploitation of the RXTE data archive is therefore especially 
important. On the other hand, sensitive radio observations 
of X-ray systems in quiescence may provide unexpected in- 
formation, and this will be easier with arrays like the Square 
Kilometre Array and its pathfinders. 
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Table 1. Fluxes, distances and outburst properties for transients 



Source v s Distance Long-term flux F av Outburst duration Outburst flux F ^ 

(Hz) (kpc) (xlO -8 ergs cm -2 s _1 ) (days) (xlO -8 ergs cm -2 s _1 ) 



Accreting millisecond pulsars 



ICR J00291+5934 


598.892 ± 2e-08 


5 ± 1 


1.8e-03 


13 


1.6e-01 


Aql X-l 


550.274 ± 9e-04 


4.55 ± 1.35 


1.2e-01 


45 


1.3 


SAX J1748.9-2021 


442.361 ± 5e-08 


8.1 ± 1.3 


9.2c-03 


61 


4.0e-01 


XTE J1751-305 


435.318 ± 4e-08 


9 ± 3 a 


2.0e-03 


10 


2.9e-01 


SAX J1808.4-3658 


400.975 ± 6e-09 


3.5 ± 0.1 


8.6e-03 


20 


3.5e-01 


HETE J1900.1-2455 


377.296 ± 5e-09 


4.7 ± 0.6 


1.8e-02 


730 


9.0e-02 


XTE J1814-338 


314.357 ± le-09 


6.7 ± 2.9 


1.3e-03 


50 


6.9e-02 


XTE J1807-294 


190.624 ± 8e-08 


8.35 ± 3.65 a 


1.4e-03 


50 


7.2e-02 


XTE J0929-314 


185.105 ± 9e-09 


7.8 ± 4.2 a 


2.7e-03 


60 


1.0e-01 


SWIFT J1756.9-2508 


182.066 ± 7e-08 


8 ± 4 a 


1.4e-04 


13 


4.0e-02 














XTE J1739-285 


1122 


7.3 ± 3.3 a 


5.9e-03 


196 


8.2e-02 


4U 1608-522 


620 


4.1 ± 0.4 


2.5e-01 


100 


2.0 


SAX J 1750.8-2900 


601 


6.79 ± 0.14 


3.5e-03 


108 


0.12 


GRS 1741.9-2853 


589 


7.2 ± 2.8 


le-02 c 


30 c 


le-01 c 


4U 1636-536 


581 


6 ± 0.5 


0.47 






X 1658-298 


567 


12 + 3 


8.0c-02 


872 


0.77 


A 1744-361 


530 b 


6 ± 3 a 


5.3e-03 


97 


0.2 


KS 1731-260 


524 


7.2 ± 1.0 


2.2c-01 


1886 


0.49 


4U 0614+09 


415 b 


2.5 ± 0.5 a 


1.2e-01 






4U 1728-34 


363 


5.2 ± 0.5 


0.23 






4U 1702-429 


329 


5.5 ± 0.2 


0.13 






MXB 1730-335 


306 b 


9.25 ± 2.85 


6.4e-02 


24 


0.59 


IGR J17191-2821 


294 b 


7.5 ± 3.5 a 


6.6e-04 


11 


0.26 


4U 1916-053 


270 b 


8.0 ± 2.2 


6.9e-02 






VD 1 oka son 




-in 1 n 

lo ± O 


o no 
y.Uc-Uz 






EXO 0748-676 


45 


7.4 ± 0.9 


4.6e-02 






cHz QPO sources 












Cyg X-2 


351 ± 34 


10.55 ± 4.45 


1.1 






GX 340+0 


343 ± 92 


9.15 ± 5.15 a 


2.8 






4U 1735-44 


316 ± 32 


8.5 ± 1.3 


8.4e-01 






GX 5-1 


288 ± 69 


7.85 ± 3.85 a 


5.3 






4U 1820-30 


285 ± 65 


7.4 ± 0.6 


1.4 






Sco X-l 


272 ± 40 


2.8 ± 0.3 


39 






GX 17+2 


272 ± 50 


11.4 ± 2 


1.8 






XTE J2123-058 


272 ± 50 


9.6 ± 1.3 


l.le-03 


53 


8.7e-02 


GX 349+2 


266 ± 13 


4 ± 0.4 


2.2 







a Distance poorly constrained: upper and/or lower limit set arbitrarily. 

b Burst oscillation frequency requires confirmation: tentative detection, or only seen in one burst from this source. 

c No reliable measurement of flux possible for this Galactic Centre source due to field crowding. The values given here are indicative 
estimates only. 



© 0000 RAS, MNRAS 000, 000-000 



Table 2. Parameter space and template requirements for the accreting millisecond pulsars. 



Source 


Us ± Aus 




P orb ± AP orb 




a x sin i ± Aa x sin i 






N Tasc 


log 10 (7Vtcmp) 


-Fstat 




(Hz) 




(hours) 




(lt-s) 




(days) 








\\jtx JUUz91+d9o4 


oyo.oyzloUoo 


18 


Z.4ob69ZZ 


23 


a. /inn'Jn no 
0.499oC-Uz 


1 


4e-07 


1 


2.6 


15.9 




± 2e-08 


(12) 


± 1.7e-06 


(15) 


± 2e-06 


(1) 




(i) 


(2.2) 


(14.7) 


Aql A-l 


£ c;n a k 

OOU.Z (45 


( .DC+UO 


1 Q C\A 7fl 

15.94/9 


1557 


2.5 


9692 


3c- 03 


1049 


15.6 


26.6 




± 9e-04 


(5.1e+05) 


± 2.0e-04 


(1038) 


± 0.5 


(6462) 




(700) 


(14.9) 


(26.0) 


bAA J174o.9-zUzl 


44z.v>blUollo 


13 


8.7b525 


139 


0.38760 


1 


4c- 06 


1 


3.8 


17.1 




± 5e-08 


(29) 


± 3c-05 


(93) 


± 4e-05 


(i) 




(i) 


(3.4) 


(15.9) 


vttt 1 TiTcri oner 


435.31799357 


31 


0.7070394 


11 


i m _ no 

l.Ulzoc-Uz 


1 


4e-7 


1 


2.5 


15.9 




± 4e-08 


(23) 


± 6e-07 


(7) 


± 5e-06 


(i) 






(2.2) 


(14.7) 


c A "V ti ono a ocrzo 

bAA J 1808.4-3658 


4UU. 975 210 2 21 


6 


z.Ulou54711 


1 


c nooi on„ no 


1 


1.0c-06 


1 


0.8 


13.2 




± 6e-09 


(4) 


± 4e-09 


(1) 


± 2.4e-07 


(i) 




(i) 


(0.6) 


(13.2) 


tlhi 1 Cj ji900.i-Z4dd 


077 Oflftl 71 f!71 

o / ( .29ol 7 19 ( 1 


5 


l.OO (0(1/ 


11 


1 o/i 1 „ no 
l.o41e-U2 


1 


7e-06 


1 


1.7 


14.7 




± 5e-09 




± 1.4e-06 


(7) 


± le-05 


(i) 




(i) 


(1.4) 


(13.2) 


Allj J 1814-338 


314.356108 /9 


9 


4.Z / 4040ZDO 


1 


0.390633 


1 


9c- 07 


1 


1.0 


13.2 




± le-08 


(7) 


± 5.6e-08 


(1) 


± 9o-06 


in 




(i) 


(0.8) 


(13.2) 


XTE J1807-294 


190.62350694 


68 


0.6678935 


1 


4.819e-03 


1 


6e-06 


1 


1.8 


14.7 




± 8e-08 


(46) 


± le-07 


(1) 


± 4e-06 


(i) 




(i) 


(1.7) 


(14.7) 


XTE J0929-314 


185.105254297 


9 


0.7263183 


1 


6.290e-03 


1 


lc-05 


1 


1.5 


14.7 




± 9e-09 


(6) 


± 8e-07 


(3) 


± 9e-06 


(i) 




(i) 


(1.2) 


(13.2) 


SWIFT J1756.9-2508 


182.065804253 


61 


0.911696 


62 


5.942e-03 


1 


6e-05 


1 


3.6 


17.1 




± 7.2e-08 


(41) 


± 2.3e-05 


(42) 


± 2.7e-05 






(i) 


(3.2) 


(15.9) 



Note. — Parameter ranges searched for the pulsars, and number of templates resulting for the mountain (r-mode) scenario assuming a coherent fold with 
T bs = 2 years. For all parameters apart from T asc we give the absolute value as well as the range searched since both quantities affect the number of templates. 
Note that AT asc is the uncertainty on the measurement, so the searched range is twice this value. 



Table 3. Parameter space and template requirements for the burst oscillation sources (details as for Table [2}. 



Source l/ 3 ±Al/ a JV„ P orb ±AP orb N P olb °x sin i ± A« x sin i N ap AT aac ^T aBC l°Slo ("temp ) -Fatal 

(Hz) (hours) (lt-s) (days) 



XTE .71739-285 



1122 

± 5 



4.2e+09 
(2.8o+09) 



120.1 
± 119.9 



2.2e+08 
(1.5c+08) 



11.3 
± 11.3 



4.1e+05 
(2.8e+05) 



1.3e+06 
(8.7c+05) 



34.4 
(33.7) 



36.0 
(36.0) 



4U 1608-522 



620 

= 5 



4.2e+09 
(2.8e+09) 



67.5 
± 57.5 



1.0e+08 
(6.7c+07) 



6.17 
± 6.17 



1.3e+05 
(8.3c+04) 



3.9e+05 
(2.6+05) 



29.1 
(28.4) 



33.7 
(33.2) 



SAX J1750. 8-2900 



601 

= 5 



4.2e+09 
(2.8o+09) 



120.1 
± 119.9 



1.2e+08 
(7.8e+07) 



11.3 
± 11.3 



2.2e+05 
(1.5c+05) 



7.0e+05 
(4.7c+05) 



33.6 
(32.9) 



36.0 
(35.5) 



GRS 1741.9-2853 



58!) 

= 5 



4.2e+09 
(2.8e+09) 



120.1 
± 119.9 



1.2e+08 
(7.7e+07) 



11.3 
± 11.3 



2.2e+05 
(1.5e+05) 



6.8e+05 
(4.6e+05) 



33.6 
(32.9) 



36.0 
(35.5) 



4U 1636-53(1 



581 

= 5 



4.2e+09 
(2.8e+09) 



3.7931263 
± 3.8e-06 



238 
(159) 



0.735 
± 8.3e-02 



1585 
(1057) 



1112 
(742 ) 



17.8 
(17.1) 



27.9 
(27.3) 



X 1658-298 



51)7 

= 5 



4.2e+09 
(2.8e+09) 



7.11610979 
± 8e-08 



3 

(2) 



1.1 
± 0.9 



1.7e+04 
(l.lc+04) 



41 

(27) 



15.4 
(14.7) 



26.0 
(26.0) 



A 1744-361 



KS 1731-2(10 



530 
= 5 



524 
± 5 



4.2e+09 
(2.8e+09) 

4.2e+09 
(2.8e+09) 



1.62 
± 0.37 



121 
± 119 



1.9e+07 
(1.3o+07) 



1.0e+08 
(6.7o+07) 



0.125 
± 0.075 



11.3 
± 11.3 



1.9e+05 
(1.3e+05) 



6817 
(4545) 

6.1e+05 
(4.1c+05) 



23.5 
(22.8) 

31.3 
(30.6) 



30.7 
(30.7 ) 

34.6 
(34.6) 



4U 0614+09 



4U 1728-34 



4U 1702-42'! 



MXB 1730-335 



IGR J17191-2821 



4U 1916-053 



XB 1254-(i!)() 



EXO 0748-676 



415 
± 5 



303 
± 5 



32!! 
= 5 



30!) 
= 5 



2!)4 
= 5 



27(1 
= 5 



95 
= 5 



45 
± 5 



4.2e+09 
(2.8e+09) 

4.2o+09 
(2.8e+09) 

4.2e+09 
(2.8e+09) 

4.2e+09 
(2.8o+09) 

4.2o+09 
(2.8e+09) 

4.2o+09 
(2.8e+09) 

4.2o+09 
(2.8e+09) 

4.2o+09 
(2.8e+09) 



0.29 
± 0.04 



5.1 
± 4.9 



120.1 
± 119.9 



120.1 
± 119.9 



120.1 
± 119.9 



0.83351411 
± 2.5e-07 



3.9334 
± 2c-04 



3.8241072 
± 2.4c-06 



3.0e+06 
(2.0e+06) 

1.5o+07 
(9.7e+06) 

6.5o+07 
(4.3e+07) 

6.0o+07 
(4.0c+07) 

5.8e+07 
(3.9o+07) 

3 

(2) 

2606 
(1737) 

9 

(6) 



7e-03 
± 7c-03 



(1.1 
± 0.1 



11.3 
± 11.3 



11.3 
± 11.3 



11.3 
± 11.3 



1.45e-02 
± 1.05e-02 



0.98 
± 0.12 



0.475 

± 0.365 



(1(1 
((,4) 

1200 
(800) 

1.2e+05 
(8.2e+04) 

l.le+05 
(7.6c+04) 

l.le+05 
(7.3c+04) 

94 
(63) 

390 
(260) 

593 
(395) 



300 
(200) 

3770 
(2513) 

3.8e+05 
(2.6e+05) 

3.6e+05 
(2.4o+05) 

3.4e+05 
(2.3c+05) 

4 

(3) 

365 
(243) 

40 
(27) 



20.2 
(19.5) 

25.5 
(24.8) 

32.9 
(32.2) 

32.8 
(32.1) 

32.7 
(32.0) 

12.2 
(11.5) 

17.7 
(17.0) 

14.5 
(13.8) 



29.0 
(28.4) 

31.7 
(31.7) 

35.5 
(35.1) 

35.5 
(35.1) 

35.5 
(35.1) 

24.0 
(23.2) 

27.9 
(27.3) 

25.3 

(25.3) 



5 



5 



3' 



("5 

s 

s 

3 



3 



CO 



Table 4. Parameter space and template requirements for the kHz QPO sources (details as for Table [2}. 



Source 










a x sin i ± Aa x sin i 


Na p 




7V TaBc 


•ogloC^temp) 


-Fstat 




(Hz) 




(hours) 




(Its) 




(days) 








Cyg X-2 


351 


2.8e+10 


236.26560000 


1 


13.5 


2.0e+04 


3.0e-02 


3228 


18.0 


27.9 




± 34 


(1.9e+10) 


± 8.3e-08 


(!) 


± 1.6 


(1.4e+04) 




(2152) 


(17.4) 


(27.3) 


GX 340+0 


343 


7.6e+10 


120.1 


8.4e+07 


11.3 


1.6e+05 


2.5 


5.0e+05 


34.5 


36.0 




± 92 


(5.1e+10) 


± 119.9 


(5.6e+07) 


± 11.3 


(l.le+05) 




(3.3c+05) 


(33.8) 


(36.0) 


4U 1735-44 


316 


2.7c+10 


4.6520042 


177 


0.693 


7825 


3.0e-3 


761 


19.1 


28.4 




± 32 


(1.8e+10) 


± 7.7e-06 


(118) 


± 0.693 


(5217) 




(508) 


(18.4) 


(27.9) 


GX 5-1 


288 


5.7c+10 


120.1 


6.9e+07 


11.3 


1.3e+05 


2.5 


4.1e+05 


34.1 


36.0 




± 69 


(3.8e+10) 


± 119.9 


(4.6e+07) 


± 11.3 


(8.7c+04) 




(2.7c+05) 


(33.4) 


(35.5) 


4U 1820-30 


285 


5.4e+10 


0.182781083 


9 


1.35c-02 


74 


2.2e-04 


28 


14.5 


26.0 




± 65 


(3.6e+10) 


± 2.8e-08 


(6) 


± 6.5e-03 


(50) 




(19) 


(13.8) 


(25.3) 


Sco X-l 


272 


3.3e+10 


18.89551 


58 


1.44 


1827 


3.0e-03 


350 


17.6 


27.9 




± 40 


(2.2e+10) 


± 2e-05 


(39) 


± 0.18 


(1218) 




(234) 


(16.9) 


(27.3) 


GX 17+2 


272 


4.2c+10 


120.1 


6.2c+07 


11.3 


1.2e+05 


2.5 


3.7e+05 


33.8 


36.0 




± 50 


(2.8e+10) 


± 119.9 


(4.2e+07) 


± 11.3 


(7.9e+04) 




(2.5e+05) 


(33.1) 


(35.5) 


XTE J2123-058 


272 


4.2e+10 


5.956833 


1657 


1.6 


3210 


1.0c-03 


422 


19.5 


28.4 




± 50 


(2.8c+10) 


± 5.6e-05 


(1105) 


± 0.3 


(2140) 




(281) 


(18.8) 


(28.4) 


GX 349+2 


266 


l.le+10 


22.5 


4.0o+05 


3.5 


3.2e+04 


0.47 


1.0e+05 


24.8 


31.7 




± 13 


(7.2e+09) 


± 0.1 


(2.7e+05) 


± 3.5 


(2.1e+04) 




(6.6e+04) 


(24.1) 


(31.2) 



Gravitational waves from accreting neutron stars 

Table 5. Computational cost effects for future performance scenarios. 



Source v X 50 improvement X 100 improvement Ta<i<ica\ 

(Hz) Attacks T ohs (days) Attacks T ohs (days) (years) 



Accreting millisecond pulsars 



IGR J00291+5934 


599 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


1.26 


Aql X-l 


550 


1 (1) 


512 


(730) 


1 


(1) 


645 


(730) 


0.17 


SAX J1748.9-2021 


442 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.35 


XTE J1751-305 


435 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.67 


SAX J1808.4-3658 


401 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.82 


HETE J1900.1-2455 


377 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.43 


XTE J1814-338 


314 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


1.12 


XTE J1807-294 


191 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.86 


XTE J0929-314 


185 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.66 


SWIFT J1756. 9-2508 


182 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


2.81 


Burst oscillation sources 




















a. i tii jT/oy-zoo 


T TZZ 


i.oe-i-uo (i.oe-rUo) 


0.1 


(0.1) 


T.Oe-f-UO 


^i.oe-i-uo ) 


U.T 


(0.1) 


n a c 

U.4o 


/ITT 1 KHQ Kilo 
4U TOUo-OZZ 


OZU 


0070 /vono'i 
00 1 z \ ( zuy ) 


3.1 


(3.6) 


/oOO 


^OODO ) 


0.0 




U.To 


civ T1 7c;n fi onnn 

o-A-A jTfOu.o-zyuu 


OUT 


i.oe-i-uo ^i.oe— |-udj 


0.1 


(0.1) 


T.Oe-f-UO 


(1 1 riA\ 
[ ( .oe-|-U4 ) 


n 1 

U.T 


j 


U.O i 


Ij-tCo T/4T.y-Zo0o 


con 


l.oe+Uo ^l.oc-l-uoj 


0.1 


(0.1) 


T.oe-f-UO 


(1 Q« 1 C\A \ 

[ ( .v>e+U4 ) 


U.T 


(U.Z ) 


U.o i 


4U 1636-536 


581 


1 (1) 


98 (167) 


1 


(i) 


123 


(211) 


0.07 


X 1658-298 


567 


1 (1) 


611 


(730) 


1 


(i) 


730 


(730) 


0.08 


A 1744-361 


530 


1 (1) 


1.2 


(2.1) 


1 


(i) 


1.5 


(2.6) 


0.61 


KS 1731-260 


524 


2.4e+04 (2.1e+04) 


0.6 


(0.7) 


2.4e+04 


(2.1e+04) 


0.6 


(0.7) 


0.08 


4U 0614+09 


415 


1 (1) 


14.9 


(25.6) 


1 


(1) 


18.8 


(32.3) 


0.31 


4U 1728-34 


363 


9.6e+04 (8.2e+04) 


3.6 


(4.2) 


9.1e+04 


(7.7c+04) 


3.8 


(4.5) 


0.11 


4U 1702-429 


329 


1.5e+05 (7.3e+04) 


0.1 


(0.2) 


7.3e+04 


(7.3e+04) 


0.2 


(0.2) 


0.14 


MXB 1730-335 


306 


7.3c+04 (7.3e+04) 


0.2 


(0.2) 


7.3e+04 


(7.3e+04) 


0.2 


(0.2) 


0.11 


IGR J17191-2821 


294 


7.3c+04 (7.3e+04) 


0.2 


(0.2) 


7.3e+04 


(7.3e+04) 


0.2 


(0.2) 


1.39 


4U 1916-053 


270 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.13 


XB 1254-690 


95 


1 (1) 


102 


(175) 


1 


(1) 


128 


(220) 


0.07 


HiJWJ U f 4o-0 / O 




1 K 1 -) 


730 


(730) 


1 

T 


(.1) 






n 17 


kHz QPO sources 




















Cyg X-2 


351 


123 (1) 


60 


(51) 


1 


(1) 


10 


(73) 


0.02 


GX 340+0 


343 


1.5e+05 (7.3e+04) 


0.1 


(0.2) 


1.5e+05 


(7.3e+04) 


0.1 


(0.2) 


0.02 


4U 1735-44 


316 


1 (1) 


35 


(61) 


1 


(1) 




(76) 


0.03 


GX 5-1 


288 


1.5e+05 (7.3e+04) 


0.1 


(0.2) 


1.5e+05 


(7.3e+04) 


0.1 


(0.2) 


0.01 


4U 1820-30 


285 


1 (1) 


730 


(730) 


1 


(1) 


730 


(730) 


0.03 


Sco X-l 


272 


1 (1) 


110 


(188) 


1 


(1) 


138 


(237) 


0.02 


GX 17+2 


272 


7.3e+04 (7.3e+04) 


0.2 


(0.2) 


7.3e+04 


(7.3e+04) 


0.2 


(0.2) 


0.02 


XTE J2123-058 


272 


1 (1) 


26 


(44) 


1 


(1) 


33 


(56) 


0.85 


GX 349+2 


266 


1908 (1625) 


41 


(48) 


1781 


(1514) 


il 


(51) 


0.05 



Note. — Computational cost effects on data analysis, assuming a 50 or 100 times improvement in computational 
power over current capabilities, for the mountain (r-mode) scenario. Where AT stac ]j. s = 1, the fold is coherent and 
T obs is the maximum feasible quantity of data that can be folded. Where Af stacks > 1, the fold is semi-coherent and 
T b a is the length of each individual data segment. The assumed spin frequency is given for each source to assist in 
cross-referencing this Table with the Figures. We also give the decoherence time T de[:o i 1 , as defined in Section [6] 
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